NVU  Nyd-9557   CI 


UNCLASSIFIED 


AEC  Computing  and  Applied  Mathematics  Center 
Institute  of  Mathematical  Sciences 
New  York  University 


TID-I4.50O  NYO-9351 

15th  Ed.  PHYSICS 

THE  SINGULARITIES  OF  THE 
RIEMANN  FUNCTION 

by 
Donald  Ludwig 

January  1,  1961 


Contract  No.  AT(30-1  }-11^80 


-  1  - 

UNCLASSIFIED 


NYO-9351 
PHYSICS  AND 
MATHEMATICS 


ABSTRACT 

This  paper  deals  with  the  Riemann  function  for 
linear  hyperbolic  systems  of  first-order  equations. 
The  leading  term  in  the  singularity  of  the  Riemann 
ftonction  is  determined  and  interpreted.   In  addition 
to  equations  with  distinct  characteristics,  certain 
equations  with  multiple  characteristics  are  treated. 
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'  •  ■  NYO-9351 

THE  SINGULARITIES  OF  THE  RIEMANN  FUNCTION 

1.   Introduction 

We  shall  deal  with  the  Cauchy  problem  for  linear  hyper- 
bolic systems  of  the  first  order,  with  n+1   independent 
variables.   Our  central  problem  is  the  description  of  the 
singularities  of  the  Rieraann  function.   In  addition  to  its 
own  intrinsic  interest,  this  problem  provides  a  key  to  the 
mathematical  theory  of  wave  propagation,  since  any  solution 
of  the  Cauchy  problem  can  be  represented  in  terms  of  the 
Riemann  function.   In  a  se-nse,  the  singularities  of  the 
Riemann  function  determine  the  structure  of  the  dependence 
of  a  solution  on  its  initial  data. 

The  analysis  of  the  singularities  of  the  Riemann 
function  also  represents  a  step  towards  an  extension  of  the 
method  of  Hadamard  to  general  linear  hyperbolic  equations. 
This  method  would  determine  the  Riemann  function  by  substi- 
tution of  a  function  of  the  proper  form  into  the  differential 
equation. 

Our  analysis  also  provides  an  approach  to  problems 
involving  propagation  of  waves  in  media  where  the  charac- 
teristics have  variable  multiplicity,  including  questions 
of  existence  and  uniqueness.   Such  problems  are  frequently 
ignored  in  the  mathematical  literature,  although  they  have 
great  physical  interest. 
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Our  principal  methods  and  concepts  are  drawn  from 
geometrical  optics  and  asymptotic  analysis.   We  shall  see 
that  the  singularity  of  the  Riemann  function  in  the  neigh- 
borhood of  a  point  is  determined  by  the  local  geometry  of 
the  ray  conoid  and  the  normal  svirface.   At  multiple  points 
of  the  normal  surface,  the  behavior  of  the  vectors  annihilated 
by  the  characteristic  matrix  must  also  be  considered. 

In  dealing  with  problems  Involving  multiple  characteris- 
tics, we  generally  restrict  ourselves  to  equations  with  constant 
coefficients  and  characteristics  with  at  most  double  multipli- 
city.—'  Our  methods  seem  to  have  an  immediate  extension  to 
more  general  problems,  but  the  complexity  of  the  method  increases 
with  the  complexity  of  the  geometry  of  the  normal  surface. 
This  is  in  contrast  to  existence  and  uniqueness  theorems 
based  on  energy  inequalities,  e.g.  for  symmetric  hyperbolic 
equations.   However,  the  same  difficulties  seem  to  appear 
when  more  detailed  Information  about  solutions  is  sought. 

We  describe  only  the  most  singular  part  of  the  Riemann 

function  at  each  point,  although  the  other  singularities  can 

2/ 

be  calculated  from  our  formulas.— '^ 

First  we  describe  our  results  for  strongly  hyperbolic 

equations  with  constant  coefficients,  whose  characteristics 

=J  We  also  must  exclude  conical  points  of  the  normal  STirface, 
or  points  where  the  characteristic  nullvectors  do  not  vary 
smoothly. 

■=/  Similar  formulas   have  b^en  obtained  by  J.  Leray 
[11],  [12],  [13]. 
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have  uniform  multiplicity.   Let  the  positive  direction  of 
the  normal  to  the  ray  cone  point  towards  the  origin.   Then 
at  regular  points  of  a  sheet  of  the  ray  cone,  the  type  of 
singularity  of  the  Riemann  function  depends  upon  the  sign 
of  the  Gauss  curvature  of  the  intersection  of  the  sheet 
with  a  plane   t  =  const.   Call  the  intersection   S. 

If  the  curvature  of   S   is  positive,  then  the  singu- 
larity of  the  Riemann  function  R   is  analogous  to  the 
singularity  for  the  wave  equation,  i.e. 

R  -x^  a  5   ^   (d) 

Here   d   is  a  measure  of  the  distance  to   S.   The  factor   a 
will  be  described  presently.   If  n   is  odd,  then  the  singu- 
larity has  support  on  S   itself.   If   n  is  even,  then  the 
singularity  has  support  on  the  positive  side  of  S. 

If  the  curvature  of  S   is  negative,  then  for   n  even 

R  -^  a  5   ^   (-d)   . 
Then  the  singularity  has  support  on  the  negative  side  of  S. 
If  the  curvat\are  of   S   is  negative  and  n   is  odd,  then 

{nil) 

R  -N^  a  log   2   Idl   , 
i.e.   R  has  a  logarithmic  singularity  on  S. 

The  factor   a   is  a  right  nullvector  of  the  character- 
istic matrix  associated  with  S.   The  magnitude  of   a   is 
proportional  to  the  square  root  of  the  Gauss  curvatiore  of 
S;   the  factor  of  proportionality  can  be  determined  solely 
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from  knowledge  of  the  tangent  plane  to   S. 

At  singular  points  of   S,   such  as  cusps,  we  are  able 
to  give  complete  results  only  if  n  =  2.   For  larger  values 
of  n,   our  methods  cover  many  interesting  cases,  and  we 
conjecture  that  the  following  description  is  correct  in  all 
cases.   Let   p-,,...,p   -,   be  the  radii  of  principal  curva- 
ture of   S.   Assume  that   p-,...,p,   vanish  at  the  point 
under  consideration,  to  orders   a-,,...,  a,   respectively. 
We  permit  any  of  the   a.   to  be  infinite.   Let 

k 


Then 


P   =  ZZ  (|  -  -^72)  . 
i=l   ^    °-l   ^ 


R  ^  a  5   "^     (d)  , 


in  the  neighborhood  of  the  singular  point  of  S. 

If  our  equation  is  weakly  hyperbolic,  i.e.  if  a  double 
root  of  the  characteristic  equation  corresponds  to  only  one 
linearly  independent  right  nullvector  of  the  characteristic 
matrix,  then  R  has  one  less  derivative  on  the  correspond- 
ing sheets  of  the  ray  cone. 

For  equations  with  variable  coefficients,  slight 
modifications  in  our  descriptions  are  required.   First  we 
consider  "planelike"  characteristic  siJirfaces   P.   These  are 
characteristics  whose  intersection  with  the  hyperplane   t  =  0 
is  an   (n-1) -dimensional  plane  passing  through  the  origin. 
The  set  of  all  such  P   envelopes  the  ray  conoid  with  vertex 
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at  the  origin.   As  before,  let   S   be  the  intersection  of  a 
sheet  of  the  ray  conoid  with  a  hyperplane   t  =  const.   At 
each  point  of  S   there  is  a  planelike  surface   P  which  is 
tangent  to   S.   Instead  of  the  principal  curvatures  of  S, 
we  consider  the  principal  relative  curvatures  of   S   and  P. 
With  this  modification,  all  of  the  results  described  above 
apply  to  equations  with  variable  coefficients.   We  also  find 
that  the  factor   a   is  proportional  to  the  square  root  of 
the  local  ray  density  on  S;   the  factor  of  proportionality 
depends  only  upon  the  bl characteristic  strip  tangent  to   S. 

Finally,  we  consider  strongly  hyperbolic  equations  with 
constant  coefficients,  whose  normal  surfaces  have  double 
points.   In  general,  a  double  point  of  the  normal  surface  is 
mapped  onto  two  rays.   We  find  that  the  Riemann  fimction  is, 
in  general,  singular  on  the  convex  hull  of  these  rays,  and 
that  the  singularity  is  given  by 

R  '^  a  6   ^   (d)  . 

Examples  to  illustrate  our  results  are  given  in 
Section  2.   For  characteristics  with  uniform  multiplicity, 
the  construction  of  the  singular  part  of  R   is  given  in 
Sections  3  and  I4..   The  method  of  stationary  phase  is 
explained  in  Section  5,  and  the  interpretation  of  results 
is  given  in  Section  6.   Equations  whose  characteristics 
have  non-uniform  multiplicity  are  discussed  in  Sections  7 
and  8. 
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This  work  has  been   stimulated  by   the   forthcoming 
Courant-Hllbert,   Methods    of  Mathematical   Physics.    II. 
I   have   benefited  from  many  discussions  with  R.    Courant, 
J.    B.    Keller,    P.    D.    Lax  emd  P.    Ungar. 

2.      Examples 

In  this  section,  we  shall  Illustrate  our  results  by 
means  of  examples.   For  simplicity,  we  restrict  ourselves 
to  equations  with  constant  coefficients. 

First  we  shall  review  some  definitions  and  ideas. 
Consider  a  linear  first-order  system  of  k   equations  with 
k  unknowns : 

.^u  s  ^A^^  +  Bu   =   0. 

Here  u  is  an  unknown  vector  function  of  the  n+1  variables 
X  ,x  ,...,x  ,  (x  =  t)  and  A  and  B  are  square  matrices 
We  assume  that   A  =1,   the  identity  matrix. 

The  characteristic  equation  is  a  restriction  on  the 
direction  of  the  normal  of  a  surface  element;  it  therefore 
refers  to  the  dual  space  with  coordinates  {X, £,..,.,,  ,£,    ), 
The  characteristic  equation  Is 

n 


det 


(XA°+  2Z  C,,  a'')   =   0  . 


v=l  ^ 


The  matrix  XA   +  >   £,      A    is  called  the  characteristic 

v=l 
matrix.   The  vectors  which  it  annihilates  are  called 

characteristic  nullvectors.   The  locus  of  points  satisfying 
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the  characteristic  equation  forms  a  cone  in  \,£,      space, 
called  the  normal  ccne  of  the  equation.   The  intersection 
of  this  cone  and  the  hyperplane   X  =  -1   is  called  the 
normal  surface .   The  ray  cone  is  the  envelope  of  all 
characteristic  planes  which  pass  through  the  origin,  i.e. 
all  planes 

"Kt  +  X  '  £,      =   0  , 

where   (X,^)   belongs  to  the  normal  cone.   The  ray  surface 
is  the  Intersection  of  the  ray  cone  and  the  hyperplane 
t  =  1. 

We  shall  be  especially  interested  in  the  relationship 
between  the  normal  surface  and  the  ray  surface.   The  simplest 
and  most  usual  situation  occurs  in  the  neighborhood  of  an 
isolated  sheet  of  the  normal  sxirface,  or  where  two  sheets 
coincide  completely.   In  this  case,  points  of  the  normal 
s\arface  are  mapped  onto  points  of  the  ray  surface  by  means 
of  the  following  construction:   Choose  a  point   P   on  the 
normal  surface,  and  draw  the  tangent  plane  "I  \     at   P.    \2^'''\f 
Let   Q  be  the  foot  of  the  normal  to  "]  ["  drawn  from  the 
origin,  and  let   S  be  the  reflection  of   Q  in  the  unit 
sphere.   S   is  the  point  on  the  ray  surface  associated 
with  P.   The  shape  of  the  ray  s"urface  at   S   is  determined 
by  the  principal  curvature  vectors  of  the  nonnal  surface  at 
P.   These  curvature  vectors  also  detennine  the  singularity 
of  the  Riemann  function  at   S, 
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normal  surface 


unit  sphere 


Figure  1 

If  P   is  a  point  where  the  normal  surface  intersects 
itself,  then  the  preceding  construction  maps   P   onto  two 
or  more  points  of  the  ray  surface.   If  the  normal  surface 
is  shaped  like  a  cone  with  vertex  at   P,   then  the  construc- 
tion maps   P   onto  infinitely  many  points  of  the  ray  siorface. 
In  general,  if  a  single  point  of  the  normal  surface  is  mapped 
onto  a  set  of  points  of  the  ray  surface,  then  the  Riemann 
function  Is  singular  on  the  convex  hull  of  that  point  set. 


A.   The  wave  equation 


Consider  the  wave  operator 


□  = 


-S 


^_i: 


Vx2 


H^      v=i  iSx")' 

The  corresponding  Riemann  function  R(t,x,tQ,XQ)   satisfies 
the  conditions: 

□  r  =   5(t-tQ,x-XQ)   =   5(t-tQ)  5(x-Xq)  , 

R  =   0  (t  <  t^)  . 
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Here   6( t-t^jX-x^),   5(t-tQ)   and   5(x-Xq)   represent  the 
Dirac  distribution  in  n+1,  1   and  n  dimensions  respectively. 
Since  the  wave  operator  has  constant  coefficients,  it  is 
sufficient  to  set   t^  ~  0*   -^0  ~  *^* 

We  recognize  that   R   is  the  solution  of  a  radiation 
problem:   the  solution  is  given  explicitly  in  [2],  p.  L(.08, 
by  setting  g(   )  =  6  (   )  .   We  shall  also  obtain  the 
solution  in  Section  l\..      Let  r  =  \x\  .      Then 


R^*^^)   =  ^"t^I-^I     """-^^  ^-°^^^ 


l|7l    ^ 


n-2  n-2      ^ 


V     r  6(t-r 

271^/^  Vdrv    *^r  vr^r;^ 


D^4.       ^      -      (-1)    ^     f     d     \^       r    6(t-r)dT;'  /  N 

R(t,x)      =     -^ n72~     — 2'  /      —  (nevenj. 


'r^-r' 

Using  the  properties   of  the    delta  function,    we  may 
combine   these  formulas    for   even  and  odd     n: 


n-1 
2 


6      2    \t2.,2) 


J   . 

■     ' 

(t  > 

0) 

k  < 

( 

(t  < 

0) 

(2.1)  R(t,x)      =    <     271 

0 

If  n   is  even,  fractional  derivatives  of  the  delta  function 

are  involved.   By  definition, 

(-1) 
e   2^/  ^  _   1    r   5(cr)do- 
o    I  s  J   -  


■/s-cr 


(s  <  0) 
(s  >  0) 


Thus  we    can  immediately  recognize   the    difference 
between  odd   and  even  dimensions;    if     n     is   odd   and     n  >   3, 
then     R     is   different  from  zero   only  on  the  ray   cone: 
t^  =  r^.      This   is   another  way  of   stating  that   the  wave 
equation  satisfies  Huygens'    principle.      If     n     is   even, 
then     R     is   different  from  zero   in  the    solid  half-cone 
t  >   r,      since   fractional   derivatives    of  the    5-function  are 
different   from   zero    on  a   half -line. 

What   we  want   to   especially  notice    about   equation   (2.1) 

Is  that      R     behaves   like    an     ^^  -  th  derivative    of  the 

2  2 

5-function  as  we   cross   the   ray    cone:      t      =  r    .      The    results 

of  Hadamard    show   that   the   Riemann  function   corresponding  to 

a   general    linear  hyperbolic   equation   of   second   order   has 

the    same    singularity   on   the   ray   conoid.      Thus   we  might    be 

led  to   conjecture   that  the   Riemann  function  always   has 

this   type    of    singularity,    if   appropriate    adjustment    is 

made    for   the    order   of   the    equation,    i.e. 

(^) 

R   ^^    5  for  first-order    systems 

(^) 
R  ^-v^-  6     "^  for    second-order  equations   or   systems 

(^) 
R  -v>    6  for    third-order   equations    or    systems, 

etc. 

Our  next  example  will  show  that  this  is  not  always  the  case, 
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- B,   The  equations  of  hydromagnetics 

Now  we  consider  the  linearized  equations  of  hydro - 
magnetics,  with  constant  coefficients.   We  shall  not  attempt 
to  construct  the  Riemann  function  for  these  equations,—^ 
but  merely  describe  the  normal  surface,  the  ray  surface,  and 
the  singularities  of  the  Riemann  function.   Our  procedures 
and  results  will  be  justified  in  Section  6. 

The  equations  for  small  hydromagnetic  disturbances 
are  as  follows:      "  


^*P0 


div  u  =   0 


Pq  J^  +   a^   grad   p    -   [i    (curl  H)   x  Hq     =      0   , 

K-r-  -   curl    (u  X  H„)      =      0 

Here   the   underlying  flow  is   at  rest,    with   constant   density 
p„,      magnetic   field   vector      H^,      and   sound  speed     a. 
-55.^  The    characteristic   equation  for    this    system   is    as 

follows : 

X[x2-(A.  C)2]    [X^-(a2+A"^)    xV   +   aV(A.  ?)^J      =      0      , 

where     A     =     /-^  H^    .      The    factor     X      is    eliminated   from  the 

— '    These   equations    are    discussed   in  Friedrichs   and  Kranzer  [5], 
H.    Grad    [6],    and  Bazer   and   Fleischman   [ij.      Recently  a 
partial    solution  for  3-dlmensional   propagation  has  been 
given  by  F.    G,    Friedlander    [i(.J,    and   a   complete    solution 
for   2-dimensional  propagation  has   been  given  by 
H.   Weitzner    [16]« 
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characteristic  equation  if  we  restrict  the  initial  data  by 

the  condition  that   dlv  H  =  0,   and  then  eliminate  one 

2    2 
unknown.   We  consider  the  case  where   A  >  a  . 

The  normal  surface  for  these  equations  is  represented 

in  Figure  2.   It  is  a  surface  of  revolution,  with  the  axis 

pointing  in  the  direction  of  Hq. 


H, 


Figure  2 

The  ray  surface  is  represented  in  Figure  3-   It  is 
also  a  surface  of  revolution.   The  outer  shell  corresponds 


H, 


Figure  3 
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to   the    inner  oval   of   the   normal    surt&ceo      The   parallel 
planes   tangent   to    the   oval   are  mapped   onto   two  points   on 
the   HQ-axis.      The    two   outer   sheets   of   the   normal    surface 
are  mapped  onto   the    cusped   triangles    inside   the    outer   shell 
of  the   ray   surface. 

First  we    consider    the    case   of    two-dimensional  propaga- 
tion,   i,eo    when  the    surfaces    are    as    shown,    and  not  revolved 
aro\md  the   HQ-axis„      Since   the   outer   shell   of  the   ray 
siorface    is   convex,    the    singularity   of   the   Riemann  function 
on  the   outer   shell    is   analogous   to   the    singularity  for  the 
wave   equation,    ioe. 

R  -^    5     '       (d)   /-^    [    0/2' 3        inside    the   ray   surface, 

R  '^    0  outside   the   ray   surface. 

Here      d     is   a  measiore    of   the    distance    to    the   ray   surface. 
The  notation      [    ^ /o-^      indicates    that   we    are    dealing  with  a 
distribution,    i.e.    that   we  must    take   finite   parts   of 
integrals    involving     R„ 

At   the    intersection  of  the    HQ-axis    and    the   outer 
shell   of   the    ray   surface,    the   Riemann  function  has   a 
sharper    singularity   due    to    the   fact    that   an  entire    sheet 
of  the   normal    surface    is  mapped   onto    a  point.      This 
situation  is   analogous   to   a   perfect   focus   in  optics.      The 
singularity   is   raised  by   a  l/2-derivative ; 

R   ^    5^^^d)    , 
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in  the  neighborhood  of  these  points.   Introducing  polar 
coordinates,  we  recognize  that  this  is  a  singularity  like 
the  delta-function  in  two  dimensions.   This  is  to  be 
expected,  since  the  ray  surface  corresponding  to  an 
ordinary  differential  equation  consists  of  a  single  ray. 
The  behavior  of  the  Riemann  function  on  the  cusped 
triangles  is  more  complicated.   The  outer  cusps  correspond 
to  the  points  of  inflection  of  the  outer  sheets  of  the 
normal  surface.   Since  the  radius  of  curvature  has  a 
non-vanishing  derivative  at  the  cusp,  the  singularity  is 
raised  by  a  l/6-th  derivative   (l/2  -  l/3  =  l/6).   This 
situation  is  analogous  to  an  Imperfect  focus  in  optics. 
Thus  we  have 

[d^/^j 
in  the  neighborhood  of  the  outer  cusps . 

On  the  curved  portion  of  the  ray  surface  between  the 
outer  cusps,  the  curvature  vector  points  away  from  the 
origin.   Here  we  must  expect  a  different  kind  of  singvilarity 
than  on  the  outer  sheet  of  the  ray  surface.   We  find: 

R  ^^    ^      '    '  {&)        on  the  side  away  from  the  origin 

Yi  '^    Q  on  the  side  toward  the  origin. 

Thus  the  singularity  is  on  the  opposite  side  of  the  surface 
than  we  would  expect  from  our  knowledge  of  the  wave  equation. 
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The  inner  cusps  of  the  ray  surface  correspond  to  the 
point  of  the  normal  surface  at  infinity.   Since  the  outer 
sheets  lie  inside  their  asymptotes,  the  ray  surface  has 
cusps o   Here  again  we  find 


R  ^^  5 


(2/3) 


(d) 


in  the  neighborhood  of  the  inner  cusps.   On  the  remainder 
of  the  ray  surface  we  have 

R  'N^  5  '        on  the  side  toward  the  origin 

R  ^^  0         on  the  side  away  from  the  origin, 

since  the  corresponding  curvature  vector  of  the  ray  surface 
points  inwardso 

We  may  sianmarlze  our  results  in  the  following  diagram: 

■1/2 


2/3^ 


Figure  \\ 
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Note  that    the   Rlemann  function   Is  not    singular   inside   the 
cusped  regions.      Weitzner's   solution   shows   that    in  fact 
the   Rieraann  function  vanishes    inside   the    cusped  regions. 
The   singularities   in  the   three-dimensional   case   may 
be   obtained  by  analogous   considerations.      We  must  bear   in 
mind,    however,    that   the   ray   s\a?face   has   an  additional 
curvature   vector  pointing  toward  the  HQ-axis.      On  the    outer 
sheet   of   the   ray   siirface,    the    singularity  is   given  by  analogy 
with  the  wave   equation: 

R  ^^    5        (d)        on  the    outer   sheet. 

However,    on  the   HQ-axis,    there   is   an  additional   point 
singularity  due   to  the   plane    sheets   of   the  normal    surface. 
At   these   points,    the   singularity  is   raised  by  one   derivative 
(1  =   2-|): 

R  ^   5^2)(^)      ^ 

in  the  neighborhood  of  these  points.   Again  introducing 
polar  coordinates,  we  recognize  the  delta-function  in  three 
dimensions,. 

The  outer  circular  cusps  of  the  ray  surface  again 
correspond  to  the  points  of  inflection  of  the  cross-section 
of  the  normal  s;ir'face<,   Here  again  the  singularity  is 
enhanced  by  a  l/6-th  derivative: 

R  ^^    5  '   (d)   near  the  cusps. 

On  the  portion  of  the  ray  surface  bounded  by  the 
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circular  cusps,  both  curvatiire  vectors  point  away  from  the 
origin.   Thus  here  again  we  have 

in  analogy  with  the  wave  equation. 

The  inner  cusps  of  the  ray  surface  correspond  to  an 
imperfect  focus  in  three  dimensions.   Accordingly,  the 
singularity  is  enhanced  by  a  l/3-derivative   (l/3  =  2(4-4)): 

R  ^^  5  ^^      id)        near  the  inner  cusps. 

On  the  remainder  of  the  ray  surface,  one  c\:irvat\are 
vector  points  inwards,  and  the  other  outwards.  Accord- 
ingly we  cannot  be  guided  by  the  wave  equation:   we  find  that 

R  -V  log^^^d)  ^    [-^J   . 

d 

This  means  that  we  must  take  finite  parts  of  integrals 
involving   l/d  .   If  we  were  to  take  logarithmic  parts  of 
these  integrals,  we  would  have  a  singularity  like   5    (d). 
Thus,  although  these  distributions  have  many  formal  simi- 
larities, they  must  be  carefully  distinguished.   In 

(2) 
particular,   log    (d)   does  not  have  support  only  at   d  =  0; 

the  Riemann  function  is  singular  inside  the  cusped  regions  as 

well  as  outside. 

We  may  summarize  our  results  in  the  following  schematic 

diagram: 
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h/l^ 


Pigiire  5 

C.   Equations  with  multiple  characteristics 

The  theory  of  hyperbolic  systems  with  multiple  character- 
istics Is  still  far  from  complete.   We  shall  confine  ourselves 
to  some  brief  remarks  and  simple  examples,  which  Illustrate 
the  theory  to  be  developed  in  Sections  7  and  8. 

If  two  or  more  roots  of  the  characteristic  equation 
coincide,  then  the  corresponding  nullvectors  of  the  character- 
istic matrix  may  not  be  linearly  Independent.   In  this  case 
the  lower  order  terms  in  the  equation  (the  matrix   B) 
determine  whether  or  not  the  Cauchy  problem  is  well  posed. 
For  example,  the  system 


(2.2) 


^t  -  ^x   =   ° 


-t   =   ° 
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is  hyperbolic,  while  the  system 

^t  -  ^x  =  0 
v^  -  u  =   0  , 

is  equivalent  to  the  equation  u. .  =  u  ,   which  Is  parabolic, 
In  Section  3,  we  shall  obtain  a  simple  condition  on  B 
which  is  necessary  and  sufficient  for  the  Cauchy  problem  to 
be  well-posed,  for  equations  whose  characteristics  have 
constant  multiplicity.   If  such   conditions  on  B  are 
necessary,  then  the  equation  is  called  weakly  hyperbolic, 
otherwise  strongly  hyperbolic.   M.  Yamaguti  and  K.  Kasahara 
[17]  have  shown  that  a  system  is  strongly  hyperbolic  if  and 
only  if  the  nullvectors  of  the  characteristic  matrices  for 
fixed  ^,,...,4   form  a  (uniformly)  complete  set. 

The  Riemann  function  for  a  weakly  hyperbolic  system  is 
more  singular  than  for  a  strongly  hyperbolic  system.   On  the 
appropriate  sheets  of  the  ray  surface,  the  singularity  of 
the  Riemann  function  is  enhanced  by  a  number  of  derivatives 
equal  to  the  defect  in  the  number  of  nullvectors  of  the 
characteristic  matrix.   Thus  the  Riemann  function  for  the 
system  (2.2)  is  given  by 

/  6{x)  t  5-^(x) 

\  0  5(x) 

In  general,  if  the  normal  surface  has  double  points, 
then  the  Riemann  function  is  singular  not  only  on  the  ray 
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STjrface,  but  also  on  certain  manifolds  which  are  supported 
but  not  envelopped  by  the  corresponding  characteristic 
surfaces.—/   Thus  the  Riemann  function  is  generally  singular 
on  the  convex  hull  of  each  sheet  of  the  ray  surface.   If  the 
normal  svirface  has  a  double  point,  then  the  singularity  on 
the  corresponding  part  of  the  hull  differs  from  what  we 
might  expect  from  consideration  of  the  wave  equation.   If 
the  equation  is  strongly  hyperbolic,  then  the  singularity 
is  weaker  by  l/2-derivative;  if  the  equation  is  weakly 
hyperbolic,  then  the  singularity  is  stronger  by  l/2-deriva- 
tive. 

For  example,  consider  the  strongly  hyperbolic  system 

u^  .  ^  .  V  =  0  , 

v^  -  V   =  0  . 


with  Cauchy  data 


u(0,x,y)   =  f(x,y) 
v(0,x,y}   =  g(x,y) 


The  solution  is 


u(t,x,y)   =  f(t+x,y)  +   7*  g(x+t-s,s+y)  ds  , 

v(t,x,y)   =  g(x,t+y)  . 

Thus  the  Riemann  function  is  given  by 

/5(t+x)  5(y)     H(-x)H(-y)  6(x+y+t) 

R  = 

V      0  5(x)   5(y+t) 

1/  This  phenomenon  is  connected  with  the  fact  that  singulari- 
ties  in  the  solutions  propagate  between,  as  well  as  alone- 
rays,  ^* 


where   H(x)   is  the  Heavislde  function  (H(x)  =  1  for 

X  >  0,   H(x)  =  0  for   X  <  0),   The  ray  surface  is  shown 

in  Figure  6o   R  has  a  singularity  like  5(d)   across  the 

segment  of  the  line   x+y+1  =  0   cut  off  by  the  x   and  y 

axeso   At  the  endpoints  of  this  segment  R  behaves  like 
the  delta-f\inction  of  two  variables,  corresponding  to  the 

fact  that  the  ray  cone  consists  of  just  two  rayso 


-1 


Figure  6 


As  an  example  of  a  weakly  hyperbolic  system,  we  consider 


u.  +  u   -  V   =   0 
t    y    X 


with  Cauchy  data 


V.  -  V   =   0  , 


u(0,x,y)   =  f(x,y) 
v(0,x,y)   =   g(x,y) 


The  solution  is 
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u  =  f(x,y-t)  +  /  g^(x,y-t+2s)  ds  , 


g(x,t+y)  , 


and  the  Riemann  function  is  given  by 


R  = 


'5(x)  5(y-t)      5^(x)  Eit^-j^) 
0  5(x)  5(t+y) 


The  ray   surface   is   shown   in  Pigiore    7«      Across   the    segment 

p  ( 1 ) 

X  =   0,      y     <   1,      R     has    a   singularity   like      5        (d),      and 

at   the   endpoints,      R     behaves   like    the   delta-function  of 

two   variables. 


5^, 


i 


Figure    7 

3,   Propagation  of  Singularities  along  Characteristics 
with  Uniform  Multiplicity 

In  this  section  we  shall  consider  the  preliminary 
problem  of  calculating  the  propagation  of  singularities 
along  regular  characteristic  sixTfaces,  where  the  normal 
cone  consists  of  isolated  sheets,  or  sheets  which  coincide 
completely.   This  problem  has  been  treated  in  D.  Ludwig  [II].] 
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for  strongly  hyperbolic  equations.   We  shall  first  give  a 
brief  summary  of  the  results  of  [114-],  and.  a  geometrical 
Interpretation  of  the  leading  term  in  the  expansion  of  the 
solution.   Next  we  shall  treat  a  special  class  of  weakly 
hyperbolic  equations  with  constant  coefficients.   In 
Section  l\.,    we  shall  show  how  the  expansions  derived  in  this 
section  can  be  used  to  construct  the  Riemann  function. 

A.   Strongly  hyperbolic  systems. 

We  consider  a  strongly  hyperbolic  system  with  variable 
coefficients: 

(3.1)     ^u  =  ^  A^x)  -^  +  B(x)  u  =   0. 


Here      u(x)      is   a   vector  with     k     components,    and     A 


V 


(v  =  0,1,..., n)  and  B  are  k  x  k  matrices  which  are 
C  functions  of  x.  We  assume  that  A  =1.  We  shall 
construct  solutions  of  (3.1)  of  the  form 

(3.2)  u  =  2Z  fJ0U))    aJ(x)  . 

j=0   J 

More  precisely,  we  wish  to  determine  the  function  ^(x) 
and  vectors   a'^(x)   (j  =  0,1,..,)   in  such  a  way  that  if 

(3.3)  ^^j^^^  "  ^j-1^^^     (j=l,2,...), 

then  (3.2)  is  a  formal  solution  of  (3«l)  for  arbitrary 
choice  of  £^{0).      Then  we  may  choose  ^o^^^      to  be  a 
distribution,  and  (3.2)  will  be  a  (formal)  weak  solution 
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of  (3.1).   since  f .{0)     becomes  as  smooth  as  we  like  for 
sufficiently  large   j,   we  may  break  off  our  series  after 
a  finite  number  of  terms  and  reduce  a  problem  involving 
singular  data  to  a  problem  with  arbitrarily  smooth  data. 
Using  the  relations  (3»3),  we  substitute  (3.2)  into 
the  differential  equation  and  collect  terms: 

Here  and  in  the  sequel,  we  shall  always  sum  over  the  index 
V   from  0   to   n.   Our  equation  will  be  formally  satisfied 
for  arbitrary  ^r)^0^      i^  w®  set  each  tertri  inside  the  brackets 
equal  to  zero  : 

(3.1+)  A^  i^^  a°   =   0  , 

(3.5)  A^  0^    J"-^    +^aJ   =   0 
From  (3.^4-),  we  conclude  that 

(3.6)  det  (A^  0^)      =  0  . 

Thus   the    surfaces     0  =   const,    must   be    characteristic.      Let 
the   rank   of     A  j^       be      k-r.      Then  there    are      r     linearly 
independent  right    and  left   nullvectors   of     A  j6    : 

^^  0^R^     =     0  p  =  l,...,r   , 

L^^  A^  ^^   =   0  Y  =!,.<.. ,r  . 

We  may  choose   L^   (y  =  1, ...,r)   such  that 
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J   u,±,..,  , 


lY.  rP      =      5^P  a,p  =   l,...,r    . 

We   also    see    from    i^-k-^    that      a       must   be    a  linear 

1  r 

rombination   of      R    ,...,R      : 

0      =     <v^      ^P    pP 


a^      =     XI     ^  ^  R 

p=l 


Setting   j  =  0   in  (3i'5)    and  multiplying  by   L'   (y  =  l,...,r), 

1  r> 

we    obtain     r      equations   to   be    satisfied   by     cT    ,  . .  . ,  cr       : 

(3.7)  L^  fZc£icr^  nh      =      0  Y  =  l,...,r. 

P=l 

We  shall  see  presently  (see  3.20)  that  this  is  a  system  of 
first  order  partial  differential  equations  with  proportional 
principal  parts.   Hence  (3.7^  can  be  solved  by  the  method  of 
characteristics. 

Once   cr  ^  has  been  determined,  we  may  obtain  a 


modulo  a  right  nullvector  of  A  0       from  (3.5) •   This 
nullvector  can  be  determined  from  a  system  of  inhomoR-eneous 

equations  analogous  to  (3.7) •   The  other  coefficients 

2   3 
a  ja-^,...,   can  be  obtained  in  a  similar  fashion. 

Now  we  can  attack  the  Cauchy  problem  for  strongly 

hyperbolic  equations.   We  impose  Cauchy  data  for  (3»l)  in 

the  form : 

(3.8)   u(0,x^,...,x^)   =  T~  f  A0{O,x^,.,.,x^))  bhx-„..,x^)    . 

j=0  J 

This  Cauchy  problem  cannot  in  general  be  solved  by 

means  of  a  single  expansion  i3»2)    since,  for  example,   b 
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may  not  be  a  nullvector  of  a  characteristic  matrix.   First 
we  construct  all  solutions  of  the  characteristic  equation 
with  initial  data 


^''{0,x^,,.,,x'^)      =  j2^(x^,...,x''} 


y   .    ,       ,  ..       .     0 

c^-   If 


If     cA     is  hyperbolic  and     x     =   t  =   0      is   spacelike,    then 
0    (x)      are    all  real    if     ^(x   5...,x    )      is   real.      If    c:^     is 
strongly  hyperbolic,    then  the   associated   right   nullvectors 
form  a   complete    set,    i.e.    any  vector  with     k      components 
can  be   represented   as   a  linear   combination  of   these   null- 
vectors.      Thus  we  write 

m       00 


(3.9)  u     =     YZYZ  ^,^0°''^^^^    a^J(x)    . 

a=l    i=0      J 


From  the  initial  condition  (3»Q),  we  see  that 

^  a''-^(0,x^,...,x'')   =  b^(x^,...,x'')  . 
a=l 

This  is  a  system  of  equations  for  the  k   scalars   cr  '^•^ 

(p  =  1,0.., r  f    o.  =   l,..,,m),   for  each  value  of   j.   Since 


a 


,ap 


the  associated  nullvectors  R  ^   (p  =  l,...,r   J  o  =  l,...,m) 
form  a  complete  set,   o~  ^"^   can  be  determined  uniquely  on 
the  initial  manifold. 


B.   Interpretation  of  the  leading  term  in  the  expansion 

By  examining  the  system  (3»7),  we  may  give  the  scalars 
cT^   a  geometrical  interpretation.   If  we  think  of  the  square 
of  cr^      as  a  measure  of  energy,  then  we  may  state:   the 
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energy  of  the  disturbance  is  transported  along  rays.-i'   More 
precisely,  cr^      is  proportional  to  the  square  root  of  the 
local  density  of  rays;  the  factor  of  proportionality  can  be 
determined  from  knowledge  of  the  first  derivatives  of  the 
phase  function  alone.   For  equations  with  constant  coeffi- 
cients, the  local  density  of  ra:'s  is  given  by  the  Gauss 
curvatvire  of  the  intersection  of  the  phase  surface  and  the 
plane   t  =  const.   This  interpretation  will  be  applied  to 
the  Riemann  function  in  Section  6. 

Our  argument  is  as  follows:   first  we  derive  the  ordinary 
differential  equation  satisfied  by  the  ray  density.   Then  a 
series  of  calculations  enables  us  to  write  the  system  (3«7) 
in  the  form 

(3.10)  ^  (<rYv7X)  ^Z:Q^^(^^vS^  =  0   Y  =  i,...,r  . 
as  P=l 

Here   s   is  a  parameter  along  the  ray,   /^\"    is  the  local 

y6 
ray  density,  and  Q'^  depends  only  upon  the  first  deriva- 
tives of  0,      i.e.   Q^^   is  a  function  defined  on  the 
characteristic  strip.   We  shall  adopt  the  summation  convention 
that  a  repeated  index   v   is  summed  from   0  to  n,   and  a 
repeated  index   i,  j  or  k   is  summed  from  1   to   n. 

In  order  to  define  the  ray  density,  we  consider  the 
canonical  system  of  ordinary  differential  equations  which 

i/This  is  a  consequence  of  the  fact  that  characteristics  have 
uniform  multiplicity.   The  solutions  given  in  Section  7  do 
not  have  this  property.   Another  such  example  is  given  in 
D.  Ludwlg  [15 J. 
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define   the  rays.      In  the   neighborhood   of   a   solution  of   the 
characteristic   equation  with  uniform  multiplicity,    we  may 
determine     jZJ,       explicitly   as   a   function  of      jZf  ,,.,., 0        ; 

(3.11)  0^     =     X(t,x,j2l^)    . 

Let     i^(t,x)      have   the    initial  values 

0{O,x)      =     ^ix)    , 

¥e  introduce  the  parameters  7  5«.»,7   in  the  initial  plane 
by  the  equations 

X   =  X(y)        i=l,,..,no 

¥e  may  choose   y  ,»..,y   in  such  a  way  that 

y""  =  H'(x)  . 

Then  we   define      x'^Csjy)      and     p^(s,y)       (v   =   O,l,o.o,n)      by 
means    of    the    system   of    ordinary   differential   equations 
(Hamilton's  equations) 


t(s,y)      =61 


(3.12) 


X  X^owoqII      « 


p.(s,y)      -     ^ f±^  i        x,..o,n 

Jx 


with  the  initial  conditions 
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t(0,y)   =   0  , 
x^(0,y)   =   X^Cy)              i  =  l,.,.,n 
P,(0,7)  =  ^M^  1  =  1 n 

Then,  since   y   is  constant  along  a  ray,  the  ratio 

given  the  ratio  of  the  ray  density  at   t  =  s   to  the  ray 
density  at   t  =  0.   Provided  that  /\  ^   0 ,      we  may  determine 
y  as  a  function  of   x,   for  fixed   t. 

In  analogy  to  the  procedure  to  obtain  the  equation 
satisfied  by  the  Wronskian  of  a  system  of  ordinary  differen- 
tial equations,  a  short  calculation  using  (3«12)  shows  that 

(3.13)  ,A   =   -  (^     l-^pn.p')  A  • 

p^x     ^i^j  dx 

Here  all  quantities  are  considered  as  functions  of  y  and   s 

Now,  thinking  of   y   as  a  function  of  x   and  t,   we 
recognize  that 

^(t,x)  =  y''(t,x)  ,  ' 

and  the  gradients  of  J   ,'»<>, J  are  tangential  to  the 

surfaces  0  =   const.   Thus  we  may  define  an  "inner  ray 
density"  as  the  (n-1 ) -dimensional  volume  spanned  by  the 
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gradients   of     J   ,  > '  -  ,7        •>      We   have 

where      p      Is   the   inner  ray  density. 

Now  we   proceed  to    analyze    the   system    (3.7)' 

(3.7)     L^^^  [A^P  i^  +   (A^ii^  +  BRP)   cr-P]     =     0. 
^  3x^  dx^ 

The  vectors  l''^   and  R^   satisfy  the  relations 

(3.15)  A^  p^  Rp(t,x,p)   =   0  ,       p  =  l,...,r  , 

(3.16)  L^(t,x,p)  A^  p^  =   0,       y=l,...,r, 

(3.17)  L^-  RP  =   5^^  . 

Recalling  that   A  =  I   and  p^  =  X(t,x,p^),   we  differenti- 
ate (3.15)  with  respect  to   p.   and  p.  : 

(3.18)  U        +  A^)  rP  +  (7.  +  A^p  J  RP   =   0  , 

Pi  -^   Pi 

(3.19)  X  ^  RP  +  (\     +A^)RP   +  U      +aJ)rP  +  (A+A^pjRp  ^=  0, 

PiPj       Pi     Pj     Pj     Pi        ^  PiPj 

Multiplying  (3. 18)  by  L^,   we  obtain 


L^  A^  rI^  =   -  6YP  X  i  =  l,...,n  ; 

hence 
(3.20)         lY  A^  rP  =   5YP  i^  V  =  0,l,...,n  . 

This  shows  that  the  direction  of  differentiation  in  (3.7)  is 
bicharacteristic.   Multiplying  (3.19)  by  '^'^0^^      aj^d  summing 
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over   i   and   J,   we  obtain 

This  equation  enables  us  to  replace  derivatives  of  R*^  by 

by  second  derivatives  of  Xo 

¥e  shall  require  one  additional  identity:   it  follows 
from  (3»11)  that 

(3o22)    0  =     -K    .    +  \        0  i  =  l,.,,,n  , 

tx^      x^  Pj   ^J 

Now  we  can  put  the  system  (3<>7)  into  the  form  (3olO)o 

We  shall  denote  partial  derivatives  of  R^  by  means  of 

subscripts  and  total  derivatives  by  means  of  the  symbol 

d/d^      }      thus 


0  ^l< 

Rf^   +  rP  0, 


It  follows  that 

I 
^^v   -  -'      -  -p^-tl  •      p.  "Ij 

+  L^  A^  RP   +  lYbrP  . 


pYP  =   lY  A^  ^  -H  lY  BR*^   =   L^rP  0^,    +  L^aJrP  0, 

^„v  p^  tl       p.  ^i: 


V 
X 

Using  (3.22)  and  (3-21),  we  have 

pYP   =   „  5^  X     0,,    +   L-^  Rf^   X  ,  +  L^  A^  R   +  L^  B  R^ 
From  (3ol3),  it  follows  that 


A     1 


pYP   =   5YP  ^  +   1    6YP  X  +  lY  RP   X    ,+lYa^^rP^  +lYb  rP 

^  p.x  Pi  X  X 
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and  hence,    using    (3»20),    (3.7)    takes   the    form 

A       '^ 
(3.23)     YZit-^^   °cr^+  J^  (T^  +  q:'^  <r^]  =  0   ,      r  =  l,.,,,r, 

p=l  ^^ 

where 

(3.214.)  Q^P  =  I  X        5yP+lYrP  X  ^  +  lVrP^+  lYbrP  . 

p.x         ^i  X        X 
Equation  (3.23)  is  clearly  equivalent  to  (3. 10). 

C.   Weakly  hyperbolic  systems 

Now  we  consider  a  weakly  hyperbolic  system  with  constant 
coeff icientSo   We  deal  with  plane  characteristic  surfaces, 
and  assume  that  only  one  characteristic  has  multiplicity  two. 
We  can  solve  the  Cauchy  problem  for  this  system  if  we  can 
exhibit  two  independent  solutions  corresponding  to  the  double 
characteristic,  since  solutions  corresponding  to  the  other 
characteristics  may  be  treated  by  the  method  outlined  for 
strongly  hyperbolic  systems. 

We  obtain  two  independent  systems  by  showing  that  the 
scalars   CT  -^   satisfy  a  second-order  differential  equation, 
instead  of  a  first-order  differential  equation.   In  our 
special  case,  the  factors  (T  ^   can  be  determined  from 
ordinary  differential  equations  along  the  rays. 

We  consider,  as  before. 


(3.25)    ^u     .     ±A^^  + 

v=0    ^x 


Bu   =   0   , 


where  A  =1.   Let   4  be  a  fixed  vector  with  components 
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^^,..,,^  .   We  can  construct  the  Riemann  function  provided 

that  we  can  obtain  solutions  of  (3»25)  which  depend  only 

0       ^   ^ 
on  X   and   >  x  £,.  .   Henceforth,  we  shall  always  sum  the 

^"^  0 

repeated  index   i   from  1   to   n.   We  denote   x   by  t, 

and   (x  ,,,.,x  )   by  x.   We  look  for  solutions  of  (3.25) 

of  the  form 

00 

(3.26)         u  =  2Z  f^(^)  aJ(t,C'x)  , 


where 


^(t,x)   =  Xt  +  C  • 


Note  that   j   is  summed  from   -1   to   en  ,   instead  of  from 
0  to   oo  .   We  assume  that  0     is  a  double  root  of  the 
characteristic  equation,  i.e.   -X   is  a  double  eigenvalue 
of  A  4 . .   We  assume  in  addition  that   A  E,.      has  only  one 
eigenvector  corresponding  to   X  ;   thus  there  exist  vectors 
R,  S,  L  and  M   such  that 

(3.27)  A^  ^^  R  +  X  R   =   0   , 

(3.28)  A^  4^  S  +  X  S   =  R   , 

(3.29)  L  A^  4^  +  X  L  =   0   , 

(3.30)  M  A^  ^^  +  X  M   =   L   . 

Substituting    (3.26)    into    {3>,2S),    we    obtain 

(3.31)  u     =      f_^{0)    [A^^a-lj    +  f.^(0)    [A^^^a°+c^-^J 

00 


+  2Z  f  ,(0)    [AXa^'""^   +c^aJ]      =      0    . 
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j=0       J  ^ 


Hence,    as   before. 


a  =     <r        R    , 


and   0-"   must  satisfy  the  equation  corresponding  to  (3.7J: 

(3.32)  L^((r"^R)   =   LA^R  cr^"-'-  +  LBR  cr"^   =   0. 

By  assijmption,  C  depends  only  upon  t   and 

X  •  ^  =  y.   Thus  (3.32)  becomes 

L  R  or""'-  +  L  A^  C.  R  cr'^   +  L  B  R  o-"-*-  =  0   , 
t  ^1    y 

But  (3.28)  and  (3-29)  imply  that   LR  =  0,   and  (3.2?)  then 
implies  that  LA^.R  =  0.   Thus  we  have 

L  B  R  o-'-'-   =   0  . 

Since  we  are  interested  in  non-trivial  solutions,  we  are 
led  to  the  condition  that 

(3.33)  L  B  R  =  0   . 

It  can  be  shown  that  (3.33)  is  a  necessary  and  sufficient 
condition  for  the  Cauchy  problem  for  (3*25)  to  be  well- 
posed.—/   If  (3*33)  is  satisfied,  then  we  can  solve  for   a 
in  terms  of   cT    : 

(AI  +  k^E,^)    a°  +  ^(cr"^  R)   =   0   , 

or 


—'     This  follows  from  A.  Lax  [lO],  and  from  an  unpublished 
note  by  A.  Lax  concerning  the  connection  between  a 
first-order  system  and  the  corresponding  single  equation 
of  higher  order. 
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(7.1 +A^4J    a°   +   (oi"-^    -   Xo-"^)   R  +  BRcr"-^      =     0   . 
1  ^  J 

Thus 

(3.31].)  -   a°   =    (cTj."^    -   X^""'")   S   +   Pcr'^  +  cr^R 

where      P     satisfies 

(3.35)  (X+A^^^)    P     =      B  R    . 

Note  that  LP  =  MBR.   Then  the  compatibility  condition 

Lci(a°)   =   0   , 
leads  to  the  equation 

(3.36)  LS(^  -  X  j^)^   cr'^   +  (MBR+LBS)  (^  "  ?^  j^)  o"  "^ 

+  LBP  o-"^   =   0  . 

Equation  (3.36)  may  be  interpreted  as  a  second-order 
ordinary  differential  equation  for  <T- 

The  coefficients   a  ,  a  ,  ...   may  be  determined  in  a 
similar  fashion  from  (3.31).   Since  (3.36)  is  a  second-order 
equation,  we  may  prescribe   cr -^   and  OT-^      independently  for 
t  =  0.   Thus  we  essentially  have  two  solutions  corresponding 
to  ^  =  \t   +   £,  '  X   ,      since   u  has  initial  values 

u(0,x)   =  f^(i2^(0,x))  [cr^R-^~^S]  +  ^  fA0{O,^))    a^  , 

if  we  prescribe  o~       (0,x)  =  0„ 

In  particular,  solutions  exist  which  have  one  less 
derivative  than  their  initial  values. 
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l\..      Construction  of  the  Rlemann  Function 

In  this  section  we  shall  show  how  the  Riemann  function 
can  be  constructed,  using  a  method  due  to  R.  Courant  and 
P.  D.  Lax  [3}'      The  singular  part  of  the  Riemann  function 
can  be  determined  by  means  of  the  expansions  derived  in 
Sections  3  and  ?•   We  shall  illustrate  the  method  by 
constructing  the  Riemann  function  for  the  wave  equation. 

We  consider  a  linear  hyperbolic  system  of  first-order 
equations : 

c^U   =   0  , 

with  Cauchy  data 

u(0,x)   =  g(x)  . 

Here   x  =  t,   and  x  =  (x  ,...,x  }. 

The  Riemann  function  R(t,x;z)   satisfies  the  conditions; 

c^R  =  0  t  >  0  , 

R(0,x;z)   =   5(x-z)  I  . 

Then  obviously 

u(t,x)   =  J*R(t,xjz)  g(z)  d(z)  . 

Thus  every  solution  of  the  Cauchy  problem  can  be  represented 
in  terms  of  the  Riemann  function.   Uniqueness  theorems  for 
hyperbolic  equations  show  that  R(t,x;z)   has  compact  support 
in  z-space  for  fixed   t   and  x. 

We  can  construct  the  Riemann  function  if  we  can 
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decompose  the  delta-f -unction  into  a  superposition  of  plane 
waves.   Such  a  decomposition  is  given  in  F.  John  [?],  p.  11: 

6(z)   = — -  J         (z*  cj)   log  (z»6J)  dcj, 

kl(27ri)^   |uj|=l 

where      log   zo6j      denotes   the    complex-valued  function  defined 
on  the    real   axis.      We    choose      k  =  1      if     n     is   odd,      k  =  0 
if     n     is   even,    and   differentiate   under   the   integral   sign. 
Thus  we    obtain 


5(z)      =      

(2 


^        /        log^^^z  -UJ)    dco    , 


where      log  denotes    the   n-th  derivative    of  the   log  function, 

in  the    sense   of    distributions. 

Now  let     w(t,xjz,6^)      satisfy   the   following   conditions: 


Xw      =      0 


t  >   0 


w 


(0,x;z,tj)      = 1_^  iog^^\(x-z)oi^)    . 


(27ri)^ 

Then 

R(t,xjz)      =         J        w(t,x;z,cj)    do;    . 
|cjI=1 

Obviously,  since  the  initial  date  are  real,  the  imaginary 
part  of   R  vanishes. 

The  most  singular  part  of   R   can  be  obtained  from  the 
expansions  given  in  Sections  3  and  8: 

(I+.l)   R  -^  ^   /   log^^^(jZ$'^(t,xiz,6J))  a"(t,x;z,6^)  dc^. 
a=l  \(J\=1 


-   kO 


Here  we  use  only  the  first  term  In  the  expansion  of  w,   in 
the  simplest  case  where  the  characteristics  are  distinct. 
Analogous  formulas  hold  for  the  other  cases.   Note  that  if 
the  operator  ^^   has  constant  coefficients,   w(t,x;z, oo) 
will  depend  only  upon  t   and   (x-z)»6J. 

The  wave  equation 

Now  we  illustrate  our  method  by  applying  it  to  the 
wave  equation.   The  Riemann  function  for  the  wave  equation 
satisfies  the  following  conditions: 


R(0,x;z)  =  0 
-4^   R(0,x;z)   =  5(x-z) 


The  corresponding  function  w(t,x;z,6j)   satisfies 
,2      n    j2 
dt^  k=l  (^x^)' 


w(0,x;z,  CJ  )   =  0 


(2TXi)^ 


^w(0,x;z,^)   =   -  -  ^  ^  log^"^((x-z).60)  . 


Since  w   and  R   are  functions  of  x-z,   we  may  set   z  =  0 
with  no  loss  of  generality.   We  shall  use  the  superscript   n 


to  indicate  that  we  have   n   space  dimensions.   Obviously 

since   w 

variables. 


since   w     satisfies  the  wave  equation  in  two  independent 
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(27ll)" 

Hence, 

2(2iii)^    lw!=l 
This   is    an  integral   of  plane  waves   over   the   unit   sphere. 

By  a  well-known  formula    [7Jj    P-    8, 

1 

(1,,3)   R^^^t.x)   =   ^  -£i^l      r  [iog(^-l^t+rp)    - 

2(2  711)'^"-^  ^1 

-   log^^"^^-t+rp)J    (l=p2)    2  dp    , 


if     n  >   2,      where 


2 


and  r  =  \x\ ,      Thus   R^^^   is  a  function  of   t   and  r 
alone.   Differentiating  our  formulaj  we  obtain 

^  =  -  :!^  /  tlog<»)(t+rp)-log'-)(-t.rp)]p(l-p2)'^dp, 

°^        2(2Tii)^  -1 

Integrating  by  parts  j,   we   have 

1  "  ' 

Ar^     _              ^  -/In-l  f    n       (n+1)/.  ,       ^    ,       (n+1)/   ^.,^„n>c 

'^j  =      -  • — — — J     [log  (t+rp}-log  (-t+rp;j^ 

^^  (n-l)2(27ii)''  "1 

n-1 

^(l-Pp)        dp    o 


Thus  we   obtain  the  well=known  formula    (see    [2],    p.    1|08 
ikA)  R^-^2)      ^      „l.iR!!^     ^ 


1^2 


5 


Thus  R^^'      can  be  obtained  for  any  n  If  we  know  R^  ' 
and  R^^K 

From  {hr.2), 

R^^^   = ^  Jr       I       [log(t+x«GJ)+log(-t+x.6j)]  dcj  . 

Introducing   0   as  an  angular  coordinate  on  the  unit  circle, 

27C 

.(2)      _     ^ 


R^'^''      =     4r  — ^   J    [loglr    cos    e+t!+log|r   cos    e-t|J    d©    . 


This    integral   can  be   found  in  tables    of   definite    integrals. 
Finally,  • 

t     >  r 


1    1           1 

r(2)      = 

^''Vt^-r^ 

0  t2  <  r2 

By  definition  of  the  fractional  Integral, 

2/i 

(  3 ) 
R  -^        is   given   immediately  by  equation   (I;.  3).      We   have 

1 

2(27ti)^  ^1 

On  the  real   axis, 

Im  log   Z      =      TC    (1    -   H(Z))    . 

Hence 

Im  log^^^    Z      =      -    Ti   5'(Z)    , 
and 

-  14-3   - 


1 

f 

2(27t)"  -1 


R^^^   =    ,'^    2  f~    -^[6'(t+rp)-5'(-t+rp)]  dp 


Thus 


R^^^   =  j^    [5(t-r)-5(t+r)]  . 
It  follows  that 

(l+,6)  R^^^   =  -^   6(t2.r2)  , 

Using  (It-olj.),  ik'S)    and  (L|..6),  we  have 
(1;.7)         R^""^   =  -^  5  ^   (t^-r^)  , 

5.   The  Method  of  Stationary  Phase 

In  this  section  we  shall  use  the  method  of  stationary 
phase  in  order  to  find  the  singularities  of  integrals  of 
distributions o   Our  results  are  formulated  in  lemmas  5.1 
and  5»2.   This  method  will  be  given  a  geometrical  interpre- 
tation in  the  following  section, 

A,   Fractional  integration  and  differentiation 
First  we  define  fractional  integration  and  differentia- 
tion of  a  distribution  which  is  singular  at  the  origin.   We 

define 

s 

(5.1)     I?  f(s)   =  — ^—  /  f(t)  (s-t)"""^  dt  , 

^  ru)  4 

A 

(5o2)     I^  f(s)   =  -~-    J   f(t)  it-s)'"'^   dt  o 

n{a)   s 


-  hh 


We  ars  only  interested  in  the  singular  part  of   I^  f,   and 
this  singular  part  is  independent  of  A   if   |s|  <  B  <  A. 
Since   s  will  always  be  restricted  to  a  finite  interval, 
we  may  formally  write  A  =  oo  ,   even  though  such  integrals 
may  not  converge.   Thus,  strictly  speaking,  we  define  the 
singular  part  of  fractional  Integrals  of   f  by  formulas 
(5«1)  and  (5»2).   In  the  sequel  we  shall  deal  only  with 
singular  parts  of  integrals. 

It  can  be  verified  in  the  usual  way  that 

ij  iP  f   =   if  P  f  , 

If  we  define  the  transformation  T  by 

T  f(s)   =  f(-s)  , 

then 

I^  f  =  T  I^  T  f  . 

I.   may  also  be  defined  for  negative   a  by  means  of 
differentiation.   Note  that 

-^  I^  f(s)   =   -  f(s)  . 
ds 

We  shall  be  especially  concerned  with  the  relationship 

1/2        1/2 
between  I.'    and  1 J     ,      Strictly  speaking,  these  operators 

do  not  commute,  but  their  singular  parts  do  commute,  i.e. 


(jl/2  il/2  .  jl/S  il/2,  f(3, 
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Is   in     C  if     f(s)      is    In     C  outside   a  closed  interval. 

Also,    taking   singular  parts   only, 

jl/2  1-1/2      ^     jl/2   i;l/2     =   y.^ 
where      /r      is   the    Hilbert   transform: 

00 

(5.3)  y   f(s)     =    /     £^^  dt  . 

-oo 

For   futia^e   reference,    we   give    some   fractional   integrals; 

a-1 


(5.1;)  ij  5(s)  =  )  r'^a) 


'   (-s)' 


s   >   0 


s   <   0   , 


s   >    0 


s   <   0    . 


If     0  <    a  <  1 


J 


{""(l-a^^cos  a  71  s*^~  s  >   0 

(5.6)      ijlog^^^(s)      = 

ni-a)    ("3)^"^  s  <  0   , 


(5.7)      I'^log^l^s)      = 


r(l-a)   s''"^  s  >   0 

-P'(l-a)cos  a7i(-s)~       3<0, 


(5o8)    (1^+1+)   log^^^s)      =      r(l-a)(l+cos  cTi)  (sgn  s)    Isl"""^, 

(5.9)  1^2    log(l^s)        =        TT    1^2     5(3)     ^ 

(5.10)  1^2    log(l-\s)        =       "TT    1^2    5(3)     . 


1^6   ~ 


B.   Two  lemmas  on  intep;rals  of  distributions 

Now  we  are  in  a  position  to  formulate  two  lemmas  regard- 
ing integrals  of  a  distribution  with  respect  to  a  parameter. 
Let   f(s)   be  a  distribution  which  is  smooth  except  at  the 
origin,  and  let   (j)(x,5)   be  a   C°°   function  of  x  and  the 
parameter  4o   We  assume  that   (4  ^,4  t**«'»4  -n^   ^  ^    '      ^^^ 

XX  X 

that  locally     4      is  determined   as   a   single-valued  function 
of     X     by  the    condition  of  stationary  phase: 

(5.11)  4^(x,4)      =      0      . 
We  write 

?     =      ^(x)    , 
and 

(5.12)  <j.(x,4)      =     4(x,^(x))   +    (C-^(x))P   iKx,C)    , 

where   \|/(x,^(x))  j^   0,   Let   +  denote  the  sign  of 
\lf(xj^(x)),   and  +  denote  its  negative.   The  integer   p 
may  depend  upon  x.   We  define  ^ix)      by 

(5.13)  4(x)   =   (|)(x,^(x))  . 

Lemma  5« 1 ?   If  a(x,4)   is  in   C   ,   with  compact  support 
in  4s   then  the  integral   J,   given  by 

J(x)  =  /  f(4(x,C))  a(x,C)  d^  , 

-00 

has  the  following  representation: 
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(a)  If      p      is   even,    then 

(5.11^)      J      =     ^  I^j+1)/P   f(^(x))   ^(x)    +  h^ix)    , 
j=0      + 

where 

aO(x)      =     ^ni^    |t(x,^)|-l/Pa(x,5)    . 

(b)  If      p      is   odd,    then 

(5.15)      J     =     y^[iiJ-^lVPf(^(x))a''>J(x)^-I^J''^^^Pf(^(x))a'4x)] 
j=0 

+     h^(x)    , 
where 

(c)  If     4(x,4)   s  $(x),      i.e.      P  =  00  ,      then   it   follows 
at   once   that 

J     =     f (J(x))    a(x)    , 

where 

a(x)      =     f  a(x,C)   dC    . 

In  any  case,   a''(x)   are  in   C°°  ,   and  the  remainder 
h^(x)   can  be  made  as  smooth  as  we  like  if   N   Is  chosen 
large  enough. 

Proof;   It  follows  from  (5«12)  that 

(|)^(x,4)   =   p(4-^)^''^[^(x,C)+(€-^)  reg.  funct.J  , 
and 
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Equation    (5ol2)    also    implies   that 


(4-^)      = 


i/P 


sgn    (4-^)    , 


and  hence 


?^x 


1  _    (3gn(g:-^))^"P        (f-I 

ITU     '  p  tix;^;       irnCTT 


III 


[l+(5-^)    reg.    funct. ] 


Now,    introducing      \      as   the    variable    of    integration,   we   have 


J      =      J^    +   J„    , 


where  ~ 

J,       =         /  (sgn(4-^))^"P   f(<^)      %   . 

1  4(x,-od)  ^   ^ 


111      P 


^        (sgn{4-^))^-P  f(c|.)   ^ 
4 


P  4* 


iz£ 
izl    P 


d(|.  , 


^ 


dcf)    . 


Here      a,(x,?)      and      ap(x,4)      are  regular  functions.      We 
note   that 


a,  (x,^) 


ap(x,^)      =     a(x,^)    . 


Now  we  assume  that   p   is  even  and  \|/(x,^)  >  0;   the 
proof  for  the  other  cases  is  obtained  by  obvious  modification. 
Since   \l/(x,^)  >  0  and   p   is  even,   <j)(x,oo)   and  (})(x,-ao) 
must  both  be  greater  than  ^.      Nothing  is  changed  if  we  set 
a-(x,4)  =   0   for   4  >  cj)(x,-oo),   and  we  set   a^lx,^)  =  0 
for   (j)  >  (|i(x,oo)<.   Thus  we  can  extend  the  limits  of  Integra- 
tion  to  the  interval   <j)  £  4'  <  oo  .   Combining  J^   and   J^  , 
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00                    ^  (1/p-l) 
(^.16)      J  =/     f(4)('l'-4)  {aAx,a   +  a.(x,4))   d(j,   . 

Now  we   expand   the    integrand   in  powers   of      <j>    : 
(5.17)     ^= tTr    (an(x,4)   +  ap(x,C))   = 

=  f-  a^(x)(M)^'^P   .  hN(,,^)    (^-$)(N-1)/P    . 

Substituting  (5.17)  into  (5»16),  we  obtain 

J(x)   =  TZ   I^J'^^^/P  f(4)  aJ(x)  +  h^(x)  , 
j=0 

which  was  to  be  proved. 

Lemma  5.1  can  be  used  to  prove  the  corresponding  result 

/I      k  \ 
when  t,      is  a  vector:   4i  =  (€  ,...,4  )•   The  most  important 

case  is  when  <j»(x)   is  a  regular  envelope  of   (j>(x,C),   i.e. 

when   p  =  2.   For  this  case,  the  more  general  result  can  be 

stated  as  follows: 

Lemma  5.2:   Let   ^  =  {^  y..,,^) ,      and  let   f(s),   (|)(x,4), 
^(x)   and  a(x,C)   satisfy  the  hypotheses  of  Lemma  5.1. 
If,  in  addition 

(j)  .  .(x,^(x))   <   0        i  =  1,...,/, 

4  .  .  (x,^(x))   >   0        l=i+l,...,k, 

(f)  .  .(x,^(x))   =   0        i  7^  j   , 
€  €^ 

then     J(x),      given  by  the    integral 
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J(x)      =     /f((j)(x,?))    a(x,.^)    dC    , 

has   the   leading  term 

li'^--^'/2lf/2f(  +  (x,5(x)))  a(x)  , 

where 

(2Tt)^^^  a(x.^(x))    ^ 

irr  4  i  ,(x,^(x))|i/2 

1=1  V£, 


•(x)   = 


This  lemma  follows  from  k  applications  of  Lemma  5ol. 

The  more  general  case,  where  certain  of  the   cj)  .  . 
vanish,  is  more  delicate.   We  shall  deal  with  such  cases 
individually. 
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6.   Geometrical  Interpretation  of  the  Method  of  Stationary 
Phase 

In  this  section  we  shall  give  a  geometrical  interpreta- 
tion to  the  method  of  stationary  phase,  as  applied  to  our 
representation  formula  for  the  Riemann  function.  We  shall 
derive  the  rules  given  in  the  introduction,  which  give  the 
singularities  of  the  Riemann  function  in  terms  of  the 
geometry  of  the  ray  conoid.   Here  we  consider  only  equations 
whose  characteristics  have  uniform  multiplicity.   The  more 
general  case  of  intersecting  sheets  of  the  normal  surface 
will  be  treated  in  Section  8. 

A.   Discussion  of  Results 

We  have  shown  in  Section  'i—'     that  singularities  in 
solutions  of  linear  hyperbolic  equations  propagate  along 
characteristics,  and  that  the  magnitude  of  the  singularity 
is  proportional  to  the  square  root  of  the  local  ray  density. 
Since  the  Riemann  function  is  a  solution  of  the  partial 
differential  equation,  we  should  expect  that  the  singulari- 
ties of  the  Riemann  function  obey  the  same  law.   We  shall 
verify  this  fact,  using  the  method  of  stationary  phase. 

We  shall  also  see  that  the  sign  of  the  determinant 
which  gives  the  local  ray  density  determines  the  character 


TT 


More   details   are   given  in  D.   Ludwig    [ll|]. 
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of  the  singularity  of  the  Rleinann  function,  since  this  sign 
determines  the  parity  of  the  integer  /O-  which  appears  in 
Lemma  5«2. 

In  the  case  of  equations  with  constant  coefficients, 
the  local  ray  density  may  be  identified  with  the  Gauss  cxorva- 
ture  of  the  intersection  of  the  sheet  of  the  ray  cone  in 
question  and  the  plane   t  =  const.   Thus  we  are  led  to  the 
rules  given  in  the  introduction,  which  relate  the  number  of 
outward  pointing  principal  curvature  vectors  and  the  charac- 
ter of  the  singularity  of  the  Riemann  function. 

Our  results  also  apply  to  equations  with  variable  coeffi- 
cients.  We  consider  "planelike"  characteristic  surfaces  which 
envelope  the  ray  conoid.   These  are  characteristic  surfaces 
which  intersect  the  initial  hyperplane   t  =  0   in  an  (n-l)- 
dimensional  plane.   Then  instead  of  the  Gauss  curvature  of 
the  ray  conoid,  we  consider  the  curvature  of  a  sheet  of  the 
ray  conoid,  relative  to  the  curvature  of  the  planelike  charac- 
teristic surface  which  is  tangent  to  the  ray  conoid  at  the 
point  in  question.   This  relative  curvature  is  proportional  to 
the  relative  ray  densities  within  the  sheet  of  the  ray  conoid 
and  the  planelike  characteristic.   With  this  modification,  our 
results  are  valid  for  equations  with  variable  coefficients: 
the  character  of  the  singularity  of  the  Riemann  function  is 
determined  by  the  number  of  inward  pointing  directions  of 
relative  curvature  along  certain  curves.   The  magnitude 
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of  the  singularity  is  given  by  the  product  of  the  relative 
curvatures  and  the  coefficient  of  the  singularity  along  the 
planelike  characteristic  surface.   Thus  the  magnitude  of 
the  singularity  is  proportional  to  the  square  root  of  the 
local  ray  density  on  the  ray  conoid. 

The  limiting  forms  of  these  expressions  near  cusps  of 
the  ray  conoid  suggest  the  rules  given  in  the  introduction 
for  determining  the  singularity  of  the  Riemann  function  at 
a  cusp.   In  its  present  form,  the  method  of  stationary 
phase  falls  to  predict  the  behavior  of  the  Riemann  function 
at  all  possible  singularities  of  the  ray  conoid^  but  those 
examples  which  are  amenable  to  treatment  by  this  method 
support  our  conjecture. 

It  should  be  noted  that  our  formulas  give  the  contri- 
bution to  the  singularity  of  the  Riemann  function  from  a 
single  point  of  stationary  phase  on  the  unit  sphere.   The 
singularity  at  a  particular  point  on  the  ray  conoid  is  the 
sum  of  the  singularities  arising  from  the  corresponding 
points  of  stationary  phase  on  the  unit  sphere.   Since  the 
point  opposite  to  a  stationary  point  is  also  a  stationary 
pointj  all  of  these  formulas  should  be  multiplied  by  a 
factor  of  at  least  two.   In  addition,  a  single  point  on 
the  ray  conoid  may  correspond  to  points  on  different  sheets 
of  the  normal  surface,  and  these  contributions  must  also 
be  added. 


Sk 


Bo   Equations  with  constant  coefficients  In  three 
independent  variables. 

First  we  shall  consider  equations  with  constant 
coefficients,  setting  n  =  2.  We  shall  discuss  only  the 
leading  term  in  the  expansion  of  the  Riemann  function, 
although  the  other  terms  can  be  obtained  in  a  similar 
fashion.   We  have  an  integral  of  the  form 

(6.1)     J  =   J       f (<|)(t,x,a))  a(t,x,a)  da  , 

|a|=l 

where   4   is  a  solution  of  the  characteristic  equation: 

4  =  A( a)t  +  a  •  X  . 

Our  representation  formula  (I4..  l)  for  the  Riemann 

(2 ) 
function  shows  that   f(4>)  =  log     \k\      ii"  the  operator 

(3) 
is  strongly  hyperbolic,  and  f(<f))  =  log    \^\      if 

4i(t,x,9)   is  a  double  root  of  the  characteristic  equation, 

with  only  one  linearly  independent  right  nullvector. 

According  to  the  method  of  stationary  phase,   J   is 

singular  only  on  the  envelope  of  the  surfaces   4=0, 

where   a  runs  over  the  unit  circle.   Let   9  be  the 

angular  coordinate  on  the  unit  circle.   We  have  a 

singularity  at  those  points   (t,x)   where 

())(t,x,eQ)   =   0  , 
4Q(t,x,9Q)   =   0  , 

-  55  - 


for  some  value  of   0^.   We  assume  first  that   c|)qq  (t,x,9_ )  /  0. 

Let   +  =  sgn  ^QQ'      Since  the  singularity  of   J  depends  on 

the  integrand  only  in  the  neighborhopd  of   (t,x,9  ),   we  may 

use  a  partition  of  unity  in  (6.1)  and  apply  Lemma  5.1,  by 

formally  allowing  the  interval  of  integration  to  be   (-ao,oo). 

Lemma  5*1  Implies  that 

-/2ia(t,x,9  )     / 
(6.2)        J  ^  ^^—   1^/2  f(<j)(t,x,9^))  . 

VUe9^t,x,9QJ|   T  0 

The  function   <})(t,x)  =  (j){t,x,9  (t,x))   is  a  solution  of 
the  characteristic  equation,  and  the  vector   a( t ,x,9^( t,x) ) 
is  a  right  nullvector  of  the  corresponding  characteristic 
matrix.   The  vector   a( t,x,9_ ( t,x) )   is  given  by  the  plane 
wave  solution  described  in  Section  i^.,  with  fixed  9,   A 
short  calculation  shows  that   I  <i'Qe("t,x,9  ( t,x) )  |  ~   is  the 
curvature  of  the  curve   (|>(t,x)  =  const,   in  the  plane 
t  =  const.,   at  the  point   (t,x). 

In  order  to  determine  the  sign  of   (j'an^tjX,  9„  ( t  ,x) ) , 
we  must  examine  the  curvatiire  vector  of  the  intersection  of 
the  ray  siorface  with  the  plane   t  =  const.   Let  us  first 
asstmie   that   o^^Qq^  *  -^0  ^  ^»   i.e.  (^.    =   ^(©q)  <  0.   Then 
if  the  curvature  vector  points  toward  the  origin  at  Xq 
(Pig.  8),  we  have 

(}»(t,x,9Q(x)  )   >   (j)(t,x,eQ(xQ) )  ,   near  Xq  . 
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<|>{t,x,eQ(x))  =  0 


cKt,x,9„(x.))  =  0 


o'^o 


Figure  8 


Hence  we  must  have   (|>QQ(t,x,  9q(x)  )  <  0,   and 

sgn  (f)QQ(t,x,eQ(x))   =   sgn  X(9q(x))  . 
By  similar  reasoning  for  the  remaining  cases  we  find: 

(6.3)      sgn  4QQ(t,x,9Q{x))   =   sgnX(9Q(x))  , 

if  the  curvature  vector  points  inwards, 

(6.14.)      sgn  <l)00(t,x,9Q(x))   =  -  sgn  X(9q(x))  , 

if   the    curvature    vector   points   outwards. 

Now  we    can  combine   our  results   and  use   Lemma  5»1« 
Let      +  =   sgn  A(9-).      -^f  the    curvature   vector  points    inwards, 

R    ^    J^_  a(t,x,9^)    1^/2   1      (2)(y)    ^ 
Using    (5.9)    and    (5-10), 


(6.5)  R  ^    -^^^  a(t,x,9    )    i;l/2   5(-H-J)    , 


or 
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y2  ^    0  • 

Here      K(t^x)      is   the    cvirvatiore    of  the    intersection  of   the 
ray   cone    and   the   plane      t   =   consto,    and  the   distribution 


f\^ 


(+(j))~^/^   vanishes  if   +0  <  0 . 

In  a  similar  fashion  we  find  that  if  the  curvature 
vector  points  outwards,  then 

(6  =  7)  R  ^    -  ^  Siit.x^Q^)    |K(t,x)'r/2  (+'J)-3/2  ^ 

If   (i)Q0(t,x,0Q)  =  0,   but  ^QQQ'.t,x,QQ)   f^   0,   then 
we  obtain  from  (5«15)  and  (5<.8)  that 

(6.8)     R^   6^  JUlZiiCyZi)       a(t..,6.)    i^i-5/3   . 

These  formulas  agree  with  our  results  for  the  wave 
equation,  and  they  imply  our  results  for  two-dimensional 
propagation  of  hydromagnetic  disturbances,,  given  in 
Section  2o 

Co   Equations  with  constant  coefficients  in  mere  than 
three  independent  variables 

We  can  extend  our  results  to  apply  to  any  number  of 
independent  variables^,  using  the  results  for  three 
independent  variables.   Our  aim  is  to  apply  the  method  of 
stationai?y  phase  to  an  integral  of  the  form 


J  = 


J       f  (4)(t,Xj,a) )  ait^x^a)  da 


|a|=l 
where   (i)(t,Xja)  =  'k{a)t   +  a  <>  Xo 
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As  before,  the  singularity  of   J   is  determined  by 
the  points  of  stationary  phase: 

(i)^(t,x,aQ)   =  0  . 

Here  (j)  denotes  the  gradient  of  <j)  with  respect  to  a, 
and  a  is  restricted  by  the  condition  that  |a|  =1.  We 
can  introduce  orthogonal  coordinates  4  ,•••»€  ^^  "the 
unit  sphere  such  that  the  matrix 

(6.9)  H  =   ((1)  .  .(t,x,a^))     (i,j  =l,...,n-l), 
is  diagonal.   Here  we  consider  <\{t,x,a{E,))      as  a  function 

of   ?o 

First  we  assume  that   H   is  non-singular.   Then 
C  ,»••,?     and 

(6.10)  <|)(t,x)  =     4(t,x;4(t,x)) 

are  determined  as  functions  of  x  ,...,x   for  fixed   t, 
and  conversely  x  ,,,.,x   are  determined  as  fionctions  of 
C-*-,  ..o,?^"''"   and   (f).   We  write   4^  =  ^^(t,x)   (  i  =  1,  .  .  . ,  n-1 ). 
Prom  the  condition  of  stationary  phase: 

(6.11)  ^4(t,xK(t,x))   ^   0        i  =  l,,..,n-l, 

we  have 

n-1       \  >^k 

(6.12)  <})  ,  .  +  m  4  ^  V  ^  =  0   (i=l,,..n-l|  j=l,...,n). 

C  x^    k=l   4^^^  dx^ 

Hence,  since   H   is  diagonal. 


-59  - 


(6.13) 


where 


c)( 


M 


)   _ 


(-1)''"^  det  (U) 
det  (H; 


(6.11+)    U  = 


,n 


4  X 


^n-l  1  •  • 


4 


^n-1  n, 


Since   (i>(t,x;4)  =  ^(4)t  +  x  •  a(^), 
that   |det(U)|  =  1.   Therefore 


it  is  easy  to  see 


(6.15) 


^(xJ 


1 
det(H) 


Since   a(^)   is  the  normal  to   4)(t,x,^)  =  const.,   it  follows 
that  (9  (^,  ^)/c)  (x)   is  lust  the  ratio  of  the  area  of  the 
spherical  image  of  a  sheet  of  the  ray  cone  to  the  correspond- 
ing area  on  the  ray  cone.   Hence, 


(6.16) 


1 
detTlT 


lK(t,x)l  |A  I   =   |K(t,x)|  , 


where   K(t,x)   is  the  Gauss  curvature  of  the  intersection  of 
a  sheet  of  the  ray  cone  and  the  plane   t  =  const. 

In  order  to  determine  the  signs  of   <!>  ^  ^  (t,xj^) 
(i  =  l,...,n-l),   we  may  use  the  results  for  n  -  2j 


since 


the  coordinate  curves  Z,     =   const.   ( i  7^  j )   are  mapped  onto 
the  lines  of  curvature  of  the  ray  cone.   We  shall  demonstrate 
this  fact  in  our  treatment  of  equations  with  variable  coeffici- 
ents.  Thus  we  find  that 
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sgn   <1>    .    ^      =      3gn  7^(a^)      if  the    corresponding 
curvature   vector   points    inwards, 

sgn  (|)    .     .      =      -   sgn  A(a^)      if   the    corresponding 
curvatiire    vector   points   outwards. 

Again,  we  let  +  =  sgn  A(aQ).  Let  m  be  the 
number  of  outward  pointing  curvature  vectors.  For 
strongly  hyperbolic   equations,     (I4.0I)    states   that 

(14.. 1)      R     ^        f       log^^M({)(t,x,a) )    a(t,x,a)    da   . 
|a|=l 

The    imaginary  part   of   this   expression  vanishes.      Applying 

Lerrana  S»2.   to   this   expression,    we   find  that    if     n     is   even 

and     m      is  even, 

n-in-2     m  n-1 

R    -  1^2  i_  2        j2   iog(n)(J)    [(2^)   2    |K(t,x)  |    a(t,x,a„)]. 
+  +  — 

Since   (n-m-2)/2   is  an  integer,  we  have 

n-2  n-1 

R  _  in-m-2  J  2  ^^1/2  log^^^  (f )  [  (271)  ^   lK(t,x)l  a(t,x,a„)]. 

—    + 

Using    (5.9)    and    (5.10),   we   have 

n-1  ^  n^l 

R    -    +   71  i^"^"2   1/      5^^-^^J')    [(271)    2      |K(t,x)I    a(t,x,aQ)], 


or 

1-n  n-1 

(6.17)      R   -    1/      5(1?)    [i""'"""^  oi(27i)    2      |K(t,x)|    a(t,x,aQ)] 

In  a  similar  fashion,  we  find  that  if  n  is  even  and 
m   is  odd. 
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1-n  n-1 

(6.18)  R-     l/      5{+'^)    [i""'^   7i(2ti)    ^      lK(t,x)l    a(t,x,aQ)]. 

If     n      is    odd,    the    real   part   of     R      is   obtained  by 
taking  the    imaginary  part   of      log        ,      since      a(t,x, a)      is 
pvire    imaginary.      Since 

log  s      =     log    |sl    +   i   Tr(l-H(s))    , 

we   have 

n-l-m  n-1 

R  --    J  I_  2        i^/^   5^^'^^J')    [(2ii)    2       |K(t,x)l    a(t,x,aQ)j. 

If     m      is   even,    this   expression   simplifies   to 
1-n  n-1 

(6.19)  R  -.    l/      5(+$)    [i'^''^~^%{2n}    ^      |K(t,x)  |    a(t,x,aQ)], 

which  is  identical  to  (6.1?).   Thus  if  m   is  even,  then  the 

singularity  of   R   is  given  by  the  same  formal  expression  for 

even  and  odd  n. 

If  m   is  odd,  then 

n-l-m  n-1 

R  -^  I  I_  2    1^2  g(n-l)(y^  ^(27i)  2   lK(t,x)I  a(t,x,aQ)J. 

Using  (5.9)  and  (5-10),  this  expression  simplifies  to 

-1-n  n-1 

(6.20)  R~  ±1+^   log  H\    [i'''"''^(27r)  ^   |K(t,x)  |  a(t,x,aQ)j< 

Notice  that   R  has  a  logarithmic  singularity  only  if  n  and 
m   are  odd.    Formulas  (6.17-6.20)  justify  the  rules  given  in 
the  introduction,  for  the  case  where  the  ray  cone  is  regular. 
In  order  to  illustrate  the  procedure  where  the  ray  cone 
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is  singular,  ioSo  where  the  normal  surface  has  a  point  of 
inflection  or  a  point  at  infinity,  we  consider  the  equations 
of  hydromagneticSo   We  have  already  discussed  two-dimensional 
propagation.   For  three-dimensional  propagation,  the  analogous 
procedure  may  be  used.   Referring  to  Figures  2  and  3  of 
Section  2,  we  see  that  on  the  circular  cusps  of  the  ray 
surface,  one  of  the  principal  curvatures  is  infinite,  while 
the  other  is  finite.   Thus  we  may  apply  Lemma  5<.l  twice; 
first  with   p  =  2,   and  then  with   P  =  3-   We  conclude  that 

(6o2l)   R  ^    Uy^+iy^)    lY^   6^^h'i)< 

—1  /p  —  1  /"^ 

x[const.  14  1  il      '^222'     a(t,x,aQ)j. 
€  4        €  €  € 

Near  the  ray  cone,   R   increases  like  \^\      '"     . 

Strictly  speaking,  our  method  does  not  apply  to  the 
singularities  on  the  inner  cusps  of  the  ray  surface  on  the 
axis  of  symmetry,  or  at  other  points  on  the  axis  of  symmetry. 
The  reason  is  that  these  points  correspond  to  isolated  multi- 
ple points  of  the  normal  surface.   Near  these  points  of  the 
normal  surface,  the  associated  right  and  left  nullvectors  of 
the  characteristic  matrices  form  a  singular  vector  field. 
In  particular,  these  vectors  have  different  limits  at  the 
multiple  points,  dewending  upon  the  direction  of  approach. 
This  is  a  degenerate  kind  of  conical  singularity  of  the 
normal  surface.   Consequently  the  vector   a(t,x  a  }   is 
singular  in  the  neighborhood  of  these  points.   Nevertheless, 
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the  limiting  forms  of  oiir  expressions  (6.19)  and  (6.20) 
suggest  the  behavior  oi   R   given  in  Section  2. 

Do   Equations  with  variable  coefficients 

We  shall  treat  equations  with  variable  coefficients  by 
means  similar  to  th^"-  ■  — d  for  equations  with  constant 
coefficients.   The  method  of  stationary  phase  can  be  applied 
as  before.   The  only  new  difficulty  arises  in  the  interpreta- 
tion of  the  numbers   4  .  .(t,x;z,a^)   (i  =  l,...,n-l),   which 

4  4 
appear  in  Lemma  5*2. 

We  recall  that,  for  strongly  hyperbolic  equations, 

equation  (I4..I)  states  that 

(li..l)   R  -^  J       log^^M(l>(t,x;z,a)  )  a(t,x;z,a)  dd  . 
la|=l 

Again,  the  singular  part  of  R   comes  from  the  points  of 

stationary  phase: 

(6.22)  <i»^(t,x;z,aQ)  =  0   , 

where  h)        denotes  the  gradient  of  ^      taken  on  the  unit 
sphere. 

First  we  shall  deal  with  the  "planelike"  characteristic 
surfaces,  given  by  the  equation 

(6.23)  (i.(t,x;z,a)   =  0   . 

Here   z   and   a  are  regarded  as  fixed  parameters.   It  will 
again  be  convenient  tc  introduce   4  ,...,4  "    as  orthogonal 
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coordinates  on  the  unit  sphere  in  the  neighborhood  of 
aQ(t,x,z).   The  function  (j)   is  a  solution  of  the  charac- 
teristic equation 


i6,2k) 


^  =   X(t,X,(|)  )   , 


with  the  initial  condition 


(6.25) 


({)(0,x;z,  a)   =   (x-z)  •  a 


As  before  in  Section  3,  we  shall  solve  (6.2I4.)  and  (6.25) 

by  means  of  ordinary  differential  equations  along  rays.   In 

the  initial  plane   t  =  0,   we  Introduce  coordinates 

y  ,..«,y   by  means  of  the  relations 

n  _ 


(6.26) 


y   =   (x-z)  •  a 
y   =   (x-z)  •  a 


n-l_   /    \ 
y    =   Ix-z )  *  a 


n-1  • 


We  introduce  the  parameter   s   by  means  of  the  equations 

t(s,y)   =  1 

*  1  /       N     _       cJX(  t  ,X.p)  . 

X  (s,yj   -   -  " \'    '^  i 


(6.27) 


Po^s^y^    =   ^  (t,x,p) 

I  Pi(s,y)   =  ^^   (t,x,p) 


^x^ 


with  the  initial  conditions 


1, . •  . ,  n 


1  =  1, . . . ,n  , 
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t(0,y)  =  0 

x^(0,y)  =  X^(y)  i  =  l,...,n 

^^•2^^     ^    P(0,y)  =  a 

Po(0»7)  =  X(0,X(y),a)  . 

Here   x  =  X(y)   is  determined  from  (6.26).   Finally,  we 


determine   y   as  a  function  of  x   and   t   from  the  equations 
(6.29) 


X   -  X  (s,yj  i  =  1, . . .,n  , 


t   =   s 
Then 

(6.30)  (|)(t,x;z,a)  =  y^(t,x)  . 

Now  we  shall  show  that 

(6.31)  ^    .(t,x;z,a)   =  y  (t,x)  1  =  l,...,n-l  , 

Equations    (6.25)    and    (6.26)    show    that    (6. 31)    is   satisfied 
for      t    =  0.      It    only  remains    to    show   that      ^    .      is    constant 
along  each  ray,    since    the    rays   are   defined  by  the   equations 
y     =   const.      Differentiating    (6.2Lt.)   with  respect   to     £,    , 


we    obtain 


n 

4      i  =        ^     ^Pj      <{>      4       J  , 

C  t  j-1        J     4  x-J 


or 


n  .  . 

(6.32)  i  .     t  +  211  <l>  1    •  x^    =    0     . 

Ct  j=l     ^  xJ 

Hence      cj)    .      is    constant   along   a  ray,    and    (6. 31)    is   proved. 


4 
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In  Section  3,    the   local  ray   density   on  the    sxirface 
4   =   const,      was   defined  by   the    ratio     d  ( <l),y  )/d  (x)  . 
Equation    (6.31)    implies   that 


(6.33) 


where 


di^     =     det    (U)      , 


(6.31;) 


u    = 


4  ]_         •  •  • 

X 

V      ^ 

4  , 

4    X 

• 
« 

*.v 

^^n-1   1    ••• 

C           X 

*?"- V  / 

Now,  following  ovir  arg\jment  for  equations  with  constant 
coefficients  we  shall  analyze  the  relationship  between  the 
planelike  surfaces  and  the  ray  conoid.   We  determine  £,   ,...,4 
as  functions  of   t,x,z   from  the  condition  of  stationary  phase: 


(6.35) 


4  . (t,x|z,^(t,x,z) ) 


0 


i  —  1, . . . ,n-l  . 


At  regular  points  of  the  ray  conoid,  this  is  certainly 
possible:  we  have 

n-1       \yk 
(6.36)       4>  i  .  -^  H  4  i  k^ 

By  a  rotation  of  the  coordinates   4  ,<•••,€    ,   we  may  insure 


=   0    i=l,...n-ljj  =],... n  . 


that  the  Hessian  matrix 
(6.37)         H  = 


'^^  1  1 
4  4-^ 


V-1,1 


V^n-1 


^e-n-l^n-l 
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is  diagonal.   At  regular  points  of  the  ray  conoid,   H 
is  non-singular.   Then  (6.36)  immediately  implies  that 

//  -Jn^        ^  (cj).^)   _   /  T  ^n-1  det(U) 

where,  as  before 

(6.39)      (i.(t,x,;z)  =  (t)(t,x;z,^(t,x,z))   . 

It  follows  immediately  from  (6.36)  and  (6.32)  that 

i  ^ 

<C  (t,x,z)   is  constant  along  rays;  hence   i(<t>»^)/^(x)   is 

proportional  to  the  local  ray  density  on  the  ray  conoid. 

Then  (6.38)  and  (6.33)  imply  that   (det(H))~-'-   is  equal  to 

the  ray  density  on  the  ray  conoid  divided  by  the  ray 

density  on  the  planelike  characteristic  which  is  tangent 

to  the  ray  conoid.   This  fact,  together  with  Lemma  5*2, 

proves  our  earlier  remark  that  the  coefficient  of  the 

singiilarity  of  R  must  be  proportional  to  the  square  root 

of  the  local  ray  density  of  the  ray  conoid.   The  factor  of 

proportionality  depends  only  upon  the  characteristic  strip 

tangent  to  the  ray  conoid. 

Now  we  shall  interpret  the  numbers  ^    .  .  ( t,x;z, ^(t,x,z ) ) 

in  terms  of  the  relative  curvatures  of  the  ray  conoid  and  the 

planelike  characteristic  siirf  ace .   Let   S  be  the  intersection 

of  a  sheet  of  the  ray  conoid  with  a  hyperplane   t  =  const., 

and  let   P  be  the  corresponding  intersection  of  the  planelike 

characteristic  surface  which  is  tangent  to  the  ray  conoid  at 

a  point  of   S.   The  equation  of   S   is 
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(6,[j.O)    ({)(t,x5z)   =   (l)(t,x;z,^(t,x,z)  )   =   0   . 

If  tQ,  Xq   is  the  point  of  tangency  in  question,  the 
equation  of   P   is 

(6ol|l)         (|>(t,x;z,^(tQ,XQ;2)  )   =  0   . 

We  introduce   (})( t,x;z, ^(t^,x^;z ) )   and  ^  (t,x;z) 
as  new  coordinates  in  x-space.   In  terms  of  these  new 
coordinates^  the  planelike  surfaces   (j)(t,x;z,4)  -   const. 
are  given  by  the  equation 

(6,14-2)    \|/(t,4,^;z,4)   =  <j)(t,x((j),4);2,4)  =  const. 

The    surface      S      is  given  by   the   equation 
(6.14.3)        4'(t,(i),^|z)      =      t(t,(l>,^;z,^).     =     0      . 

It   follows   from    (6.35)    and    (6.14.2)    that 
(6.104.)        ^   1  (tp^j^jz,^)      =      0  i   =  l,..,,n-l    . 

Hence,    differentiating  with  respect   to      ^"^    : 
(6.I4.5)        ^    .     .(t,(l>,^jz,^)    +  \|/  (t,(i),^;z,^)      =     0 


i , j   -  1 , . . . , n-1 


It  follows  from  (6<,l4_3)  and  (6.l4.[|.)  that 


(6,i+.6)   \|/  ,  (t,(t),^;z)   =  \|f  .  (t,<l),^,-z,^)     i  =  l,...,n-l  . 
Differentiating  (60^.6),  we  obtain 
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1, j  =  1,. ..,n-l 


or ,  us  ing  ( 6 .  I4.5 ) , 


(6.[i.7)  ?.  .(t,(t),^jz)   =  t  i  {t,^,K}z,K)-^   i  .(t,(l),^;z,^) 

r^^  r^^  €^^ 

It  follows  immediately  from  (6.I4.2)  that  ^    .     .   =   ^    .     .. 
Hence  the  (diagonal)  Hessian  matrix  given  in  (6.3?)  repre- 
sents the  relative  ciorvatures  of   S   and   P,  measured  along 
the  (n-1)  coordinate  lines  K,     =   const.   ( i  7^  j ) .   If  the 
differential  operator  has  constant  coefficients,  then  the 
cvirvature  of   P   is  zero,  and   ^  ,...,^     from  an  ortho- 
gonal system.   Equation  (6,I|.7)  states  in  this  case  that  oior 
coordinate  lines  are  the  lines  of  principal  curvature  of  S. 

Now,  applying  the  method  of  stationary  phase  as  we  did 
for  equations  with  constant  coefficients,  we  obtain  (6.1?- 
6.20),  where   +  =  sgn  (j)^  ( t,x;z,  ^)   and   |K(t,x)|  =|det(H)r 
Singular  points  on  the  ray  conoid  may  be  treated  in  a 
similar  fashion. 
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7o   Uniform  Expansions  near  Intersections  of  the  Mormal 
Surface  o 

This  section  will  be  concerned  with  preliminary  calcu- 
lations, directed  towards  a  construction  of  the  Rlemann 
function  for  equations  which  have  self-intersections  of  the 
normal  surface.   When  two  roots  of  the  characteristic  equa- 
tion approach  each  other,  then  in  general  the  expansion  of 
the  singular  part  of  the  solution  changes  character,  i.e. 
a  "Stokes  phenomenon"  appears. 

¥e  shall  attack  the  simplest  problems  of  this  sort  by 
a  method  which  seems  to  have  wide  applicability.—   The 
method  is  based  upon  the  belief  that  the  form  of  the 
asymptotic  expansion  of  the  solution  of  a  differential 
equation  depends  only  upon  the  geometry  of  the  relevant 
characteristics o   It  follows  that  once  we  have  an 
asymptotic  expansion  of  the  solution  of  a  (simple)  problem, 
then  the  form  of  this  expansion  will  be  valid  for  all 
problems  with  a  similar  geometry  of  characteristics. 

In  the  present  case,  consideration  of  a  simple  proto- 
type leads  to  expansions  with  a  typical  term  of  the  form 


^  2  ,1 


0 


f   f(4i  +  i-^t-)  g(t,x;f)  d-e  . 


This  method  is  frequently  used  in  dealing  with  asymptotic 
expansionso   See,  for  example,  J.  B.  Keller  [8]  and 
R„  Eo  Langer  [9Jo 
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-1  p 

Where      4     7^  <|>    ,      this   expression  is    analogous    to   the 
expansions    of  Section  3,   with  wave   forms      f,  (({>    )      and 

f,  ((|)    ).      As   before,      f,       is   the    indefinite    integral   of 

1  P 

f .      However,    where      ^      =   cf)        the   wave   form   is  given  by 

A.      A   simple   example 
We   consider   the   equation 

or  the    equivalent   first-order   system: 


(7.2) 


^Jt-^i^-     =     - 


(4  -  4)  V  =  0 


^Jt  -  37 


Here  we  assume  that   X   depends  upon  certain  parameters; 

in  particular  X  =  1   for  certain  values  of  these  parameters 

We  prescribe  Cauchy  data  of  the  form 

u(0,x)   =   0 


(7.3) 


v(0,x)   =  2Z  i'-(x)  b^(x)  > 
j=0     J 


where ,    as   before,      -r-  f.(x)    =  f.   ,(x)      (i   =   1,2,...). 
'  *      dx      J  j-1  J  f    f 

From   Duhamel ' s    principle   we   find  that 

00 
v(t,x)      =     2Z  ^^(t+x)    a^(t,x)    , 

and 
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where 


t 

u(t,x)  =  J    v!{t,x;zr)   dv    , 
0 


{4r-  -  \  4-)   w(t,x;^'  )  =  0 


w 


( 'c'  ,x;  t)     =  v(  r',x)  . 


Hence  w  has  the  form 

OP  ,. 

w(t,x;r)  =  2Z  fv(^(t-r)  +'C+x)  gJ(t,x;^)  , 

and 

t 

00    n  i 

(7oi4.)  u(t,x)  =  ^  J  f  .(Xt+x+(l-X)r)  g'J(t, x;tr)  dr. 

j=0  0  J 

This  expansion  for   u  exhibits  the  behavior  described  In 

the  introduction  to  this  section.   We  have  only  to  observe 

that  two  solutions  of  the  characteristic  equation  are 

(|)-'-(t,x)   =   Xt  +  X   , 

<t^(t,x)   =   t  +  X    , 
and 

1.  X  =  i-i   . 


B.   Strongly  hyperbolic  systems. 

Now  we  consider  a  general  first-order  system  with 
constant  coefficients: 

(7.5)     <^u  =  u.  +  ^  A^  ^  +  Bu  =   0   . 
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We  permit  two  roots  of  the  characteristic  equation  to 
coincide  for  certain  directions  of  the  wave  normal,  but 
require  that  the  corresponding  characteristic  nullvectors 
vary  smoothly  and  remain  linearly  independent. 

Guided  by  the  solution  ij  .l\.) ,    we  look  for  solutions 
of  (7.5)  of  the  form 

00         J  00        . 

(7.6)  u(t,x;a)   =  7~  f  (<}>)  a''(t,x;a)  +  y~  t A^}hHt.xia) 

j=0  J  j=0  J 

00  t 

+  2Z  /  f.(<|)+P^)  gJ(t-r,x;a,r)  d-cf 
j=0  0   J 
OP  t^ 

+  211  /  ^^(^-ps)  h^(t-s,x;a,3)  ds   . 
j=0  0  J 

Here   <j)(t,x;a)   and  \l/{t,x;a)   are  "plane  wave"  solutions 
of  the  characteristic  equation: 

(7.7)  <l>(t,x;a)   =  \^{a)t  +   a»x      , 


(7.8)         \l/(t,x;a)   =  >L^(a)t  +  a 


and 


•  X   , 


(7.9)  p(a)   =  ^^(a)  -  X^ia)      . 

We  make  the  simplifying  assumption  that  the  quantities 
^    f    ^^>  g"^*  h"^   do  not  depend  upon  x   or   a.   The  construc- 
tion of  Section  l\.   shows  that  the  Riemann  function  can  be 
represented  in  terms  of  such  solutions  of  {7.$}. 

Let 
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(7.10)  a'^*     =     X^(a)   I   +  ZI  A^    a,      , 

1=1  ^ 

(7.11)  A*     =     ^^(a)    I   +  XI  A^   a,      ; 

i=l  ^ 

thus     A*''      and     A        are    the    characteristic  matrices    corres- 
ponding to     <j)      and     ^f      respectively.      Ovr   assumptions 
about    (^     iinply   that   vectors      R^,    R^,    L^      and     L^      exist, 
which  depend   smoothly  on      a,      such  that 

A  V      =      0    ,         Lh^      =      0    ,         L^  •  R^      =      1    , 
(7.12)<aV      =      0    ,         lV      =     0    ,         L^.  R*      =      1    , 
L^.R^   =      0    ,         L^.R'^   =      0      . 

Now  we   apply    (7*5)   to    (7.6),    and   collect   terms: 

(7ol3)      0  =  ^u  =   f_^((^)[A^a°j  +  f_^(^)    [aVj 

+  ^  f  .((^)    [kK^'^^+Xiah    +  h^(0;t)] 
j=0      J 

+  ^f.(f)     [A^bJ'''^+<^(b"^)    +   gJ(0;t)J 
j=0     J 

+    y  f  ,(<j)+p'Z^)   A*^   g°(t--r-'C)    d-cr 
0      "^ 

+    Jf  -.(il'-ps)   A^  h°(t-s;s)    ds 
0      '^ 

j=0  t)  J 

+  ^/  [A^h>5''^(t-s;s)+c;ifhJ(t-s;s)]   f,(^-ps)   ds    . 
j=0   0  J 

Thus,  in  analogy  to  the  procedure  of  Section  3,  we  are  led 
to  the  conditions: 
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(7.1i|)  A*^  a°(t)  =   0   , 

(7.15)  A^  b°(t)  =   0   , 

(7.16)  A<^a-^"'^(t)  +  /a^(t)  +  h^(0;t)  =  0   (j  =  0,l,...), 

(7.17)  A*bJ''^(t)  +  ^bJ(t)  +  g^(0;t)   =  0   (j=0,l,...), 

(7.18)  A^  g^{t-V}t)  =   0   , 

(7.19)  A*  h°(t-r';ir)  =  0   , 
t 

(7.20)  /f,(<t.+p'rf  [A'^gJ''^(t-r;r)  +  ^ght-Vft-)]   dr 

n   J 


t 

+  /  f  (\|f-ps)[A*hJ''^-^(t-3;s)+A^(t-s;s)]  ds  =  0  , 
0  J 

(j  =  0,1,...). 

It  should  be  noted  that  the  essential  difficulty  (which 
we  shall  presently  overcome)  is  that  the  inverses  of   A^   and 
A^   do  not  depend  continuously  on   a  ;   when   p(a)   is  small 
it  is  not  sufficient  to  require  that   L*''z   =  0   or   L^z  =  0 
in  order  to  have  a  bounded  solution  to 

A^  y  =   z   , 
or 

A*  7  =   .   , 

respectively.      However,    if  we   require   that      L^z   =  L   z   =  0, 
then  both  equations   have  bounded   solutions. 
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Over   equations    (7.16)    and    (7.1?)    thus    can  be    solved  if 

we  require   that 

L^LXa^    +  hh0}t)]      =      0      , 


(7.21) 


(7.22) 


L^[^a^    +  hJ'(0;t)]      =     0      , 
L^EXb^    +  g^(0;t)]      =      0      , 

L^[.€bJ   +   gj'(0;t)J      =      0      . 


The    integrals   in    (7.20)    can  be    treated   in  a    similar  fashion 

if  we   set      s  +  V   =   t    ;      then  we   obtain 

t 

(7.23)   /  f.(4+P'c)    [A*g^"'^(t-r;tr)   +Xg^(t-r;r) 
0     J 

+  A^h^'^-^it:'  }t-x^)    +  dChh^  ',t-^)]    dX       =      0    . 

It  should  be  noted  that  the  operator  S/dt  acts  on  the 
first  argument  of  hMt:^  ;t- 'C  ) .  Thus  we  are  led  to  the 
condition 

(7.21+)  k^s^'^^it-T^i'!:^)   +-^gJ(t--r;-r)  +  A'''hJ"^^(  ^;  t- r^ ) 

+  o^h^(r';t-'r )  =  0  . 
There  is  no  loss  of  generality  if  we  set 

(7.25)  A*hJ'''^(^;t-^)   =   0   . 

In  order  to  solve  (7.2L|.)  for   g''   ,   we  impose  the 
conditions 
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(7.26) 

L^  [o6ght-T^}^)   +  ^hU^it-^)]     =     0   . 

Now  we  may  solve  (7.114.-7.19),  (7.21),  (7.22)  (7.21+- 
7.26)  in  a  recursive  fashion.   We  shall  show  how  the  first 
terms  in  (7.6)  are  obtained;  the  others  may  be  obtained  in 
a  similar  manner.   Prom  (7.114-)  and   (7.15),  we  have 

(7.27)  a°  =  0--L  R*^  , 

(7.28)  b°  =  cr^R^  ; 
similarly  (7.18)  and  (7.19)  imply 

(7.29)  g°  =  ^^l  R^  , 

(7.30)  h°  =  y^2  ^"^   ' 

Using  (7.12),  we  may  determine   o^   and  cr^      from.  (7.21) 
and  (7.22): 

(7.31)  L*^  c2f(cr^R'^)   =   0   , 

(7.32)  L^  o^(c^R^   =   0   . 

Here  the  initial  values  for   cr^   and  cr^     may  be  prescribed 
arbitrarily.   Once   cr,   and  cr^      are  known,  we  may  determine 
ti^(0;t)   and   ^2^°**^   ^^°"^  (7-21)  and  (7-22): 
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(7.33)        L*  [^a°  +  M.2(0;t)R^]   =   0   , 

Then  a   and  b   may  be  determined  within  a  multiple  of 
R*^   and  R^   respectively  from  (7.16)  and  (7.17). 

Now  we  can  attack  the  system  (7.26).   We  introduce 
new  variables  £,      and  h       by  means  of  the  relations 

4  =  t  -  -r-  , 
Using  (7.12),  il .2b)   becomes 

^^  +  l'^'br^i  +  i}m^^^    =  0  , 


(7.35) 

L^BR^^X3_  ^^^^^   L^BrV^   =   0   . 

The  lines   C  -   const.   and   ^  =  const.   are  the  charac- 
teristics of  this  system.   |i.p   is  prescribed  for   ^  =  0 
by  (7.33).,  and  \i^      is  prescribed  for  E,   =  0      by  (7.31^). 
Thus   IX,   and   [Xp  may  be  obtained  by  solving  a  charac- 
teristic Initial  value  problem  involving  only  two 
independent  variables.   If  we  did  not  assume  that   a%  b", 
g%   and  h*'   are  independent  of  x   and   a,   then  the 
ni:imber  of  variables  involved  would  be  larger. 

Once  (7.35)  has  been  satisfied,  we  can  determine 
g-'-(t-r';'C^)   modulo  R^   from  (7.214.).   Continuing  in  a 
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similar  fashion,  we  may  obtain  a%  b%  g*^   and   h-^   in 
succession  if  initial  data  for   a^   and  b*^   are  piven. 
For   j  >  0,   it  will  be  necessary  to  solve  inhomogeneous 
ordinary  differential  equations  analogous  to  (7.3I)  and 
(7.32)  and  inhomogeneous  characteristic  initial  value 
problems  analogous  to  (7.35),  (7.33),  (7.3I4.). 

The  entire  expansion  for   u   is  determined  uniquely 
as  soon  as  Cauchy  data  for   u   are  given.   The  assignment 
of  initial  values  for   a*^   and  b*^   can  be  made  exactly 
as  in  Section  3. 
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8o   The  Slnpijularity  of  the  Riemann  Function  on  the 
Hull  of  the  Ray  Cone  . 

In  this  section,  we  shall  use  an  expansion  of  the 
Riemann  function  which  is  valid  even  for  equations  with 
self -intersections  of  the  normal  surface.   We  can  apply 
the  method  of  stationary  phase  to  integrals  of  the 
expansions  derived  in  the  previous  section.   The  dis- 
tinctive feature  of  the  Riemann  function  in  the  present 
case  is  that  the  singularity  of  the  Riemann  function  is 
not  confined  to  the  ray  cone.   Indeed,  if  a  single  point 
of  the  normal  surface  is  mapped  onto  two  rays,  then  the 
Riemann  function  is,  in  general,  singular  on  the  convex 
hull  of  these  rays.   For  strongly  hyperbolic  equations, 
the  singularity  on  the  hull  is  weaker  by  l/2-derivative 
than  it  is  at  ordinary  points  of  the  ray  cone. 

First  we  observe  that  the  decomposition  of  the 
delta-ftinction  given  in  Section  I4.,  together  with  the 
results  of  the  previous  section,  enables  us  to  represent 
the  singular  part  of  the  Riemann  function  as  a  super- 
position of  plane  waves.   In  order  to  obtain  the  highest 
singularity  of  the  Riemann  function,  we  may  apply  the 
method  of  stationary  phase  as  described  in  Sections  5  and 
60   The  only  new  feature  arises  in  the  treatment  of 
integrals  of  the  form 

(8,1)    J        J   f((j)(t,x;a)  +  p(a)T:')  git-x:',^)    dr  da  , 
|al=l   0 
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where   p(a)   may  vanish  on  certain  curves  on  the  unit 
sphere. 

The  condition  of  stationary  phase  means  that,  at  the 
point   Cq,  r-Q   (0  <  TT-Q  <  t). 


(8.2) 


4^(t,x;ao)  +  T-Q  p^ia^)     =  o  , 


(8.3)  p(aQ)   =  0   . 

Here  the  subscript   a  denotes  the  gradient  on  the  sphere 
Ia|  =  1.   Qq   and  ^^      depend  upon  t   and  x.   Equation 
(6.3)  shows  that,  at  the  point  of  stationary  phase, 
({)(t,x;aQ)  =  \|/(t,x;aQ),   i.e.   a^   corresponds  to  a  double 
point  of  the  normal  surface. 

In  order  to  interpret  T^q,      we  first  find  x,   and 
Xp   such  that 


c  f 


^^{t,x^,a^)      =     0  ,      (j.(t,x^,aQ)   ^ 

t^(t,X2,aQ)   =   0  ,      iKt,X2,aQ)   =   c   . 

The  points  x^   and  x^   are  uniquely  determined  by  these 
equations.   If   c  =  0,   x^   and  x^      represent  the  two  rays 
which  correspond  to   a^.   Since   (f)   and   t   are  linear  in 
X,   the  unique  solution  of  the  equations 

<|)(t,x,aQ)   =  c   , 

<t>^(t,x,aQ)  +  ^oP^(aQ)   =   0  , 
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is  given  by 


0        0 


Thus  T"^  is  a  parameter  on  the  segment  joining  x-,  and 
Xp.  If  X  does  not  lie  on  this  segment,  then  there  will 
be  no  point  of  stationary  phase. 

For  simplicity,  we  shall  treat  only  the  most  usual 
case,  where  the  sheets  of  the  ray  cone  corresponding  to 


<f)(t,x;a)   and   \|/{t,x;a)   are  convex  at   x^   and  Xp 

respectively.   Then,  if  £,,..,,£,  are  local  coordinates 

1        ? 
on  the  unit  sphere,  and   X  ( a^ )  =  X  ( a„ )  >  0,   the  matrices 


(f)  .  .(t,x,,a(4))  i,j  =  l,...,n-l 


and 


^    .    ,(t,X5,a(4))  l,j  =  l,...,n-l 

are  positive  definite.   Then  it  follows  that  if 

X   -   (1  -  — ^  ^1  ■"  T  ^2   * 
then 


It  follows  that  the  matrix 


^.i.j^^'^'^0^  ■*■  '^oP.i.J^^O^      ^'J  "  l,...,n-l 
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]        n— 1 

is  positive  definite.   We  may  choose   4  ,...,4     such  that 
this  matrix  is  diagonal. 

Now  we  may  apply  the  method  of  stationary  phase  to 
(8.1),  integrating  first  with  respect  to   a.   It  is  easily 
seen  that 

(8.5)  -^  [(l>(t,x;a(4(^))  +  r  p(a(c(r}))] 

P  1 


n-l 
1 


i-1  ^^x^x   P^i^x 


Hence,  applying  Lemma  S-2.   and  setting  _+  =  sgn  4^  > 
2-n 


(806) 


5(+(t)) 


i^"2  it(27t)^/2  g(t-rQ;rQ) 


n-l 


1=1 


*^i^l*^0P^i^l 


1/2 


n 


(7 


i^  '*  i.i''V^i.i 


y2 


Here  ^  =   (j>(t,x;a„).   The  same  formula  holds  for  both  odd  and 
even  n.   It  should  be  remembered  that  this  term  arises  only 
if  X  lies  on  the  segment  joining  x,   and   Xp,   i.e.   only 
on  the  hull  of  the  rays  corresponding  to   a„. 
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